Best 

Available 

Copy 


AD-783  605 


EXTENDED  ARRAY  EVALUATION  PROGRAM. 
SPECIAL  REPORT  NO.  11.  EVALUATION  OF  THE 
NORWEGIAN  SHORT-PERIOD  ARRAY 


Frode  Ringdal,  et  al 

Texas  Instruments,  Incorporated 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dm  Entered) 


03 


REPORT  DOCUMENTATION  PAGE 


1  REPORT  NUMBER 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


2.  GOVT  ACCESSION  NO.  3.  RECIP’CnT'S  CATALOG  NUMBER 


Ti  fL  E  (and  Subtitle) 


S  TYPE  OF  REPORT  A  PERIOD  COVEREO 


EVALUATION  OF  THE  NORWEGIAN  SHORT- 
PERIOD  ARRAY  -  FINAL  REPORT 


Final 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7  AuTHORfA, 


ALEX(01  )-STR-73-l  1 


B  CONTRACT  OR  GRANT  NUMBERfO 


Frode  Ringdal  and  Richard  L.  Whitelaw 


F336<5f  -72-C-0725 


PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 


Texas  Instruments  Incorporated 

Equipment  Group 

Dallas,  Texas  7S222 


10  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  ft  WORK  UNIT  NUMBERS 


VELA  T/2705/B/ ASD 


II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 


12  REPORT  DATE 


Advanced  Research  Projects  Agency 
Nuclear  Monitoring  Research  Office 
Arlington,  Virginia  22202 


M  MONITORING  AGENCY  NAME  ft  ADDRESSfn  dltlerent  tr om  Controlling  Ottice) 


2  November  1973 


13  NUMBER  OF  PAGES 

94 


IS.  SECURITY  CLASS,  (o t  thlt  report) 


Air  Force  Tecnnical  Applications  Center 
VELA  Seismologic al  Center 
Alexandria,  Virginia  22  314 


16.  DISTRIBUTION  STATEMENT  (ot  thlt  Hepott ) 


UNCLASSIFIED 


1 5a  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


APPROVED  FOR  PUBLIC  RELEASE,  DISTRIBUTION  UNLIMITED 


17.  DISTRIBUTION  STATEMENT  (ol  th»  abstract  entered  In  Block  20,  II  dlllerent  Irom  Keport) 


’  NATION  At  TFCHNIOAt 
INFORMATION  SERVICE 

U  t>  [>p.itt'nent  of  Commerce 
Springfield  VA  ?r!M 


iB  supplementary  notes 


ARPA  Order  No.  1714 


IS  KEY  WORDS  (Continue  on  reverie  aide  il  necea  sary  and  Identity  by  block  number) 


Seismology  Signal-to- noise  improvement 

NORSAR  short-period  array  Detection  threshold 
Signal  processing  Discrimination  capability 

Array  beamforming  Maximum  -  likelihood  estimation 


20  ABSTRACT  (Continue  on  reverse  side  II  •  ecestary  and  identity  by  block  number) 


This  report  describes  the  final  phase  of  tne  evaluation  of  the 
Norwegian  short-period  Seismic  Array  (NORSAR)  which  has  been  conducted 
by  Texas  Instruments  Incorporated  at  the  Seismic  Data  Analysis  Center. 

The  report  also  summarizes  results  achieved  during  the  complete  evaljation 
study  over  the  period  1  April  1971  to  30  September  1973. 

The  major  areas  of  study  covered  by  this  report  are: 


DD 


1473  EOITION  or  1  NOV  6S  IS  OBSOLETE  , 


‘  _ UNCLASSIFIED 

I  1  SECURITY  CLASSIFICATION  OF  THIS  ►  AGE  (When  Dele  Fnte-ed) 


J 


UNCLASSIFIED 


aECUWiTV  CLASSIFICATION  OT  THU  P  AQei’H?t«n  Dtlm  Ent *r»d) 


20.  continued 


•  Array  processing  performance 

•  Partial  array  performance 

•  Maximum  likelihood  estimation  of  the  NORSAR  event 

detection  capability 

•  Evaluation  of  short-period  discriminants. 

The  accumulated  data  base  for  this  study  has  been  547  events, 
including  30  presumed  explosions;  33  of  which  are  from  the  Eurasian 
continent. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAOEfltTi»n  D ala  Enfred) 


APPROVFD  FOR  PUBLIC  RELEASE,  DISTRIBUTION  UNLIMITED 


\ 


ALEX(01)  STR  73  11 


EVALUATION  OF  THE  NORWEGIAN  SHORT-PERIOD  ARRAY  -  FINAL  REPORT 


SPECIAL  REPORT  NO.  11 

EXTENDED  ARRAY  EVALUATION  PROGRAM 


Prepared  by 

Frode  Ringdal  and  Richard  L.  WhUelaw 


TEXAS  INSTRUMENTS  INCORPORATED 
Equipment  Group 
Post  Office  Box  6015 
Dallas,  Texas  75222 


Prepared  for 

AIR  FORCE  TECHNICAL  APPLICATIONS  CENTER 
AFTAC  Project  No.  VELA  T/2705/B/ASD 
Alexandria,  Virginia  22314 


Sponsored  by 

ADVANCED  RESEARCH  PROJECTS  AGENCY 
Nuclear  Monitoring  Research  Office 
ARPA  Program  Code  No.  2F10 
ARPA  Order  No.  1714 


2  November  1972 


Acknowledgment:  This  research  was  supported  by  the  Advanced 

Research  Projects  Agency,  Nuclear  Monitoring  Research  Office 
under  Project  VELA-UNIFORM,  and  accomplished  under  the  tech¬ 
nical  direction  of  the  Air  Force  Technical  Applications  Center 
under  Contract  No.  F33657-72-C-0725. 


Equipment  Group 


ABSTRACT 


This  report  describes  the  final  phase  of  the  evaluation  of  the 
Norwegian  short-period  Seismic  Array  (NORSAR)  which  has  been  conducted 
by  Texas  Instruments,  Incorporated  at  the  Seismic  Data  Analysis  Center,  The 
report  also  summarizes  results  achieved  during  the  complete  evaluation  study 
over  the  period  1  April  1971  to  30  Stptember  1973. 

The  major  areas  of  study  covered  by  this  report  are: 

•  Array  processing  performance 

•  Partial  array  performance 

•  Maximum  likelihood  estimation  of  the  NORSAR  event  detection 
capability 

•  Evaluation  of  short -period  discriminants. 

The  accumulated  data  base  for  this  study  has  been  567  events, 
including  39  presumed  explosions;  33  of  which  are  from  the  Eurasian  continent. 


Neither  the  Advanced  Research  Projects  Agency  nor  the  Air  Force 
Technical  Applications  Center  will  be  responsible  for  information  contained 
herein  which  has  been  supplied  by  other  organizations  or  contractors,  and 
this  document  is  subject  to  later  revision  as  may  be  necessary.  The  views 
and  conclusions  presented  are  those  of  the  authors  and  should  not  be  inter¬ 
preted  as  necessarily  representing  the  official  policies,  either  expressed  or 
implied,  of  the  Advanced  Research  Projects  Agency,  the  Air  Force  Technical 
Applications  Center,  or  the  US  Government. 
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SECTION  I 


INTRODUCTION 


This  report  presents  the  results  of  the  final  phase  of  the  eval¬ 
uation  of  the  short  period  Norwegian  Seismic  Array  (NORSAR).  It  extends 
the  analysis  reported  in  Special  Report  No.  6,  (Barnard  and  Whitelaw,  1972) 
and  Special  Report  No.  9,  (Ringdal  and  Whitelaw,  1973)  under  the  Extended 
Array  Evaluation  Program.  The  overall  objectives  of  the  NORSAR  SP  eval¬ 
uation  have  been: 

•  Determine  the  best  processing  methods  for  enhancing  the  signal 
to-noise  ratio  of  Eurasian  events 

•  Determine  the  array  detection  capability  for  Eurasian  events 

•  Evaluate  the  performance  of  short  period  discriminants  at 
NORSAR 

•  In  conjunction  with  long  period  NORSAR  data,  determine  the 
detection  and  discrimination  capability  of  NORSAR  for  Eurasian 
events. 

The  fourth  objective  stated  above  will  be  the  topic  of  a  forth¬ 
coming  report  (Special  Report  No.  13  ,  1973),  and  will  not  be  discussed  here. 
Five  analysis  tasks  were  undertaken  in  order  to  meet  the  first  three  objectives 

•  Noise  analysis 

•  Signal  analysis 

•  Evaluation  of  array  processing  effectiveness 

•  Detection  threshold  estimation 
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Analysis  of  the  behavior  of  SP  discriminants. 

Final  results  from  noise  and  signal  analyses  were  presented  in 
Special  Reports  No.  6  (1972)  and  No.  9  (1973).  No  additional  work  on  these 
subjects  has  been  undertaken  since  then.  Our  efforts  during  the  final  phase 
have  been  focused  on  increasing  the  data  base  for  detection  and  discrimination 
analysis  and  develop  methods  for  reliable  estimation  of  the  array  capabilities. 

Section  II  of  this  report  describes  the  results  from  the  array 
processing  evaluation,  including  a  study  of  the  performance  of  a  potentially 
reduced  NORSAR  array.  Section  111  presents  regional  estimates  of  the  NORSAR 
detection  thresholds,  while  SP  discrimination  capabilities  are  discussed  in 
Section  IV.  An  extensive  summary  of  ail  major  results  achieved  during  the  SP 
evaluation  program  is  presented  in  Section  V. 

The  NORSAR  SP  array  is  centered  about  100  km  due  north  of 
Oslo,  Norway,  at  a  latitude  of  60.  8°N  and  a  longitude  of  10.  8°E.  The  array 
consists  of  132  short  period  seismometers  and  has  an  aperture  of  about  100  km. 
The  sensors  are  grouped  in  22  six-element  subarrays;  each  subarray  has  a 
center  sensor  and  a  five  sensor  ring  and  is  7-1C  km  in  diameter  (Figure  1-1). 

The  results  presented  in  this  report  are  based  on  seismic  events 
and  presumed  explosions  from  1971  and  1972.  A  total  of  567  events  have  been 
processed;  all  but  15  of  these  were  from  E  irasia.  The  number  of  presumed 
explosions  totals  33  from  Eurasia  and  6  from  the  Western  Hemisphere.  The 
complete  data  base  for  the  NORSAR  SP  evaluation  is  listed  in  Appendix  A;  jt 
includes  the  344  events  analyzed  in  Special  Reports  Nc.  6  (1972)  and  No.  9 
(1973)  as  v/ell  as  an  additional  223  events,  mostly  from  June  and  July  of  1972, 
which  have  been  processed  since  then. 

Figure  1-2  presents  a  breakdown  of  the  processed  events  by 
information  source.  Our  events  have  been  selected  from  seismic  bulletins 
provided  by  four  different  organizations: 


1-2 


a  01C  (9) 
a  D 

a  n  n  a02C  (10) 


a  n 

n  Q  14C  (22) 


□ 

□  _  □ 


03C  (11) 


°  D  13C  (21) 

o  n 


01 B  (2)  □ 

n  □ 


n  a  07 B  (8) 

D  ° 


□  D  n 

n  02B  (3) 

□  n  '  ' 


u  03B  (4) 

a  n 


□  u  04 C  (12) 

□  □ 


a  n 
u  □  n 


n  05 C  (13) 


□  a  J2C  (20) 


a  n  06B  7) 
a  '  ' 
o 


a  D  □01A(1) 


n  a 

„  u  p06c<14> 


p 

□  n 

a 

p  11C  (19) 

n 


10C  (18) 

a  ' 


04B  (5) 


n  a 


D  a  05 B  (6) 


„  09C  (17) 

n  a 


□  a  a 

□  n08C(16) 


a  07 C  (15) 

n 


FIGURE  1-1 


20  km 


NORSAR  SHORT  PERIOD  ARRAY 


wmmmvA 

»y////i 


S1N3A3  3D  339UF1N 


[UWVI 


w 

< 

10 

os 

< 

o 

os 

2 

Q 

< 

0 

a 

2 

co 
1— 1 

>s 

>s 

>> 

>s 

JO 

o 

JO 

JO 

T) 

TO 

TO 

TO 

<b 

OJ 

0) 

0) 

+•> 

*-» 

+-» 

♦-* 

«-l 

Ll 

L» 

Li 

0 

0 

0 

0 

a 

CL 

CL 

a 

<u 

ai 

<u 

<U 

os 

os 

OS 

os 

1-4 


•  The  PDE  listings  (Preliminary  Determination  of  Epicenters) 
issued  by  the  ’‘'ational  Oceanic  and  Atmospheric  Administra¬ 
tion  (NOAA). 

•  The  LASA  seismic  bulletin  issued  by  the  Seismic  Data  Analysis 
Center  (SDAC). 

•  The  NORSAR  seismic  bulletin  compiled  at  Kje Tier,  Norway. 

•  The  bulletin  from  the  International  Seismic  Month  (ISM),  which 
covers  February  20  to  March  19,  1972,  and  has  been  compiled 
at  Lincoln  Laboratories. 

In  all  cases  where  an  event  was  reported  by  more  than  one 
organization,  we  selected  our  source  information  according  to  the  priority 
list:  ISM,  PDE,  LASA,  NORSAR.  As  can  be  seen  from  Figure  1-2,  most 
events  of  m^  above  4.  5  have  been  reported  by  PDE,  while  LASA  and  NORSAR 
supply  the  detection  information  for  most  low  magnitude  events. 

Our  procedure  in  analyzing  these  events  was  first  to  estimate 
the  P-wa''e  signal  arrival  time  at  NORSAR  for  each  event,  and  then  to  re¬ 
quest  SP  data  tapes  from  the  NORSAR  Data  Center  covering  intervals  of  ap¬ 
proximately  10  minutes  around  each  expected  arrival.  These  tapes  were 
then  processed  using  our  NORSAR  Short  Period  Array  Evaluation  Software 
Package  (Texas  Instruments  Incorporated,  1971)  to  determine  signal  para¬ 
meters  and  to  produce  plots  of  waveform  traces.  Finally,  for  each  event, 
the  TI  analyst  made  a  decision  as  to  whether  or  not  the  event  could  be  detect¬ 
ed  on  NORSAR  data,  based  upon  visual  inspection  of  the  signal  tracea  and 
personal  judgement. 

Figure  1-3  is  a  histogram  showing  the  distribution  of  process¬ 
ed  events  as  a  function  of  epicentral  distance  from  NORSAR.  The  data  base 
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may  conveniently  be  split  into  three  regions,  where  the  majority  of  events 
are  located; 


•  The  Mediterranean  area  (distance  15-30  degrees).  Approx¬ 
imately  100  events  from  this  region  have  been  processed. 

•  Iran  and  Central  Asia  (distance  35-50  degrees).  The  number 
of  events  is  approximately  150. 

•  Kamchatka  and  the  Kurile  Islands  (distance  60-70  degrees). 

Approximately  200  events  from  this  region  are  included  in  our 
data  base. 

As  can  also  be  seen  from  Figure  1-3,  very  few  of  the  process¬ 
ed  events  were  not  detected  on  NORSAR  SP  data,  with  the  exception  of  the 
Kuriles-Kamchatka  area  events. 


Throughout  1971  and  1972,  the  quality  of  the  NORSAR  SP  data 
was  excellent.  For  a  total  o.  25  of  the  events  selected  by  TI,  no  data  was 
available  from  NORSAR.  This  corresponds  to  less  than  five  percent  of  our 
data  requests,  and  thus  indicates  that  the  array  was  operational  for  an  aver¬ 
age  of  more  than  95  percent  of  the  time  For  about  one-third  of  the  process¬ 
ed  events  all  132  sensors  were  operational.  In  most  other  cases  one  or  two 
subarrays  were  dead  or  contained  calibration  signals;  the  worst  data  loss  was 
33  sensors.  Data  spikes  were  observed  for  ten  events,  but  in  each  case  only 

one  or  very  few  sensors  were  affected;  consequently  these  events  could  still 
be  satisfactorily  processed. 


The  phase  reversals  observed  in  1971  data  (Special  Report 

No.  9,  197  3)  appeared  to  have  been  permanently  corrected  from  about 
January  1972. 


SECTION  II 


ARRAY  PROCESSING  PERFORMANCE 

A.  INTRODUCTION 

In  the  course  of  the  NORSAR  SP  evaluation  study,  several  array 
processing  techniques  were  examined  for  possible  application  to  NORSAR  short 
period  data.  It  was  decided  at  ar,  early  stage  that  adaptive  multichannel  methods 
would  not  be  very  effective,  since  the  seismic  noise  is  essentially  incoherent  in 
the  frequency  band  on  interest.  Although  signal  frequency  characteristics  show 
large  variations,  we  found  that  a  "standard"  filter  (Special  Report  No.  6,  >972) 
would  in  most  cases  give  close  to  optimum  signal-to-noise  ratio  (SNR)  improve¬ 
ments.  Also,  significant  time  delay  anomalies  (deviation  from  plane  wave  pro¬ 
pagation  along  the  great  circle  azimuth)  were  observed  for  all  regions,  thus 
causing  cons:derable  beamforming  loss  when  plane  wave  dela-  .,  were  applied 
between  subarrays.  Finally,  we  observed  large  signal  amplitude  variations 
from  subarray  to  subarray;  this  suggested  that  some  form  of  weighted  beam¬ 
forming  might  be  beneficial. 

We  thus  concluded  that  the  following  array  processing  techniques 
should  be  evaluated  in  detail  on  NGRSAR  data: 

•  Bandpass  filtering,  using  the  standard  filter 

•  PJar.e  wave  subarray  beamforming 

•  Adjusted-delay  array  beamforming 

•  Diversity-stack  array  beamforming  (i.  e.  ,  applying  weights  pro¬ 

portional  to  subarray  signal  amplitudes  prior  to  beamforming). 
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ihis  section  presents  the  results  from  these  analyses.  The 
data  base  contains  statistics  from  412  Eurasian  events,  of  which  172  occurred 
in  1971  and  the  remainder  in  the  first  six  months  of  1972. 

Also  included  in  this  section  is  a  brief  discussion  of  the  per¬ 
formance  of  a  potentially  reduced  NORSAR  array.  The  most  natural  way  to 
reduce  the  size  of  the  array  is  to  eliminate  the  outer  ring  of  subarrays;  thus 
leaving  8  subarrays  as  shown  in  Figure  II- 1.  However,  the  Northeastern 
^.art  of  the  cr^y  has  consistently  shown  the  best  signal-to-noise  ratios  for 
Eurasian  events;  therefore,  we  decided  also  to  evaluate  the  performance  of  a 
partial  array  consisting  of  8  subarrays  in  the  Northeast  corner  as  shown  in 
Figure  11-2.  A  total  of  60  Eurasian  events  from  1971  were  selected  as  data 
base  for  the  comparison  of  these  tvo  partial  arrays  to  the  f"!l  array. 

In  the  subsequent  discussions,  beamforming  and  filtering  gains 
have  been  obtained  by  computing  the  signal-to-noise  ratio  (SNR)  for  each  event 
at  each  processing  stage,  and  then  defining  gain  as  SNR  improvement  from  one 
stage  to  another.  The  following  formula  has  been  applied: 

SNR  =  10  •  log,„  ( 

&10  N  ' 

where  S  is  the  signal  plus  noise  power  over  a  6.  4  seconds  window;  (this 
window  length  is  comparable  to  signal  duration  for  medium  size  events  and 
is  also  convenient  for  computational  purposes;  representing  a  power  of  two 
relativ<  to  our  sampling  rate  of  0.  1  seconds).  Similarly,  N  represents  the 
noise  power  averaged  over  a  window  of  approximately  one  minute  preceding 
the  onset  of  the  signal. 

B.  ARRAY  BEAMFORMING  AND  FILTERING  PERFORMANCE 

NORSAR  short  period  beamforming  is  implemented  as  a  two 
step  process.  The  first  step  is  subarray  beamforming,  wnich  is  performed 
separately  for  each  of  22  subarrays.  The  second  step,  which  we  will  refer 
to  as  array  beamfurming,  consists  of  adding  together  22  cubarray  beams  with 
the  proper  time  delays. 
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FIGURE  II- 1 

GEOGRAPHICAL  LAYOUT  OF  A  PARTiAL  NORSAR  ARRAY 
(FILLED  CIRCLES)  CONSISTING  OF  THE  A  AND  B  RINGS 
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FIGURE  II-2 

GEOGRAPHICAL  LAYOUT  OF  A  PARTIAL  NORSAR  ARRAY 
(FILLED  CIRCLES)  SITUATED  IN  THE  NORTHEAST 
CORNER  OF  THE  FULL  ARRAY 
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Tho  pe rformance  of  NORSAJ1  subarray  beamforming  was  eval- 
uated  in  Special  Report  No,  9,  (1973)  and  it  has  been  verified  that  plane  wave 
delays  in  general  are  appropr  ate  for  this  purpose.  SNR  gains  for  subarray 
beamforming  averages  approximately  7  dB. 

Table  II - 1  shows  means  and  standard  deviations  of  array  beam 
forming  and  standard  filter  SNR  gains  for  those  events  of  1971  whose  average 
subarray  SNR  is  at  least  5  dB.  The  statistics  from  this  set  of  events  are  con 
sidered  representative  of  ideal  array  pe  'formance,  since  the  beams,  in  this 
case,  were  generated  using  adjusted  deJayj  determined  by  examination  of  the 
subarray  bearn-to-reference  subarray  beam  c  ros  s -cor  relation  functions.  In 
contrast,  the  array  beam  of  each  1972  event  was  generated  using  delays  de¬ 
termined  from  the  previously  processed  master  event  with  the  nearest  epi¬ 
center.  Statistics  from  these  events  can  be  regarded  to  determine  a  conser¬ 
vative  estimate  of  operational  array  performance. 

The  following  principal  points  are  noted  in  regard  to  Table  11-1 
and  the  additional  statistics  discussed  above: 

•  Array  beamforming  gain  from  average  subarray  beam  to 
adjusted-delay  beam  is  1  0  dB  +  2  dB  for  both  wide-band  and 
filtered  data  (Figure  II-3).  This  compares  with  a  theoretical 
(VN)  value  of  13  dB.  Outliers  on  the  low  end  of  the  distribu¬ 
tion  shown  in  Figure  11-3  represent  events  from  epicentral  re¬ 
gions  such  as  Turkey,  Greece,  and  Western  Russia,  which 
typically  suffer  high  beamforming  losses  as  a  result  of  poor 
signal  similarity  across  the  array. 

•  Reference  subarray  to  adjusted-delay  beam  SNR  gain  is  not 
very  large.  Events  from  1971  averaged  4.  5  +  3  dB  for  wide¬ 
band  data  and  slightly  less  for  filtered  data.  The  distribution 
of  these  gains  is  shown  in  Figure  11-4.  Outliers  on  the  lew 
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TABLE  II- 1 


NORSAR  .SHORT  PERIOD  ARRAY  PERFORMANCE  STATISTICS 
FOR  172  EVENTS  FROM  1971 


Beamforming  Gains  ] 

Wide-Band 

Standard  Filter 

Average  Subarray  to 

Adjusted-Delay  Beam 

10.  2  ±  2.  1  dB 

9.  8  +2.  1  dB 

Reference  Subarray  to 

Adjusted-Delay  Beam 

4.  6  +  3.  1  dB 

3.  6  +  2.  8  dB 

Adjusted-Delay  Beam  to 

Diversity-Stack  Beam 

1.  0  +  0.  9  dB 

1.  6  ±  1.  0  dB 

|  Standard  Filter  Gains 

Reference  Subarray  Beam 

8.6  +4.3  dB 

Average  Subarray  Beam 

8.  1  +  3.  7  dB 

Adjusted -Delay  Beam 

7.  8  +4.  2  dB 

Diversity-Stack  Beam 

8.  3  +  4.  1  dE 

Combined  Gain 

(Wide-Band  Average  Subarray  Beam  to  Filtered  Adj 

usted-Delay  Beam) 

17.  9  ±  4.  5  dB 
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DISTRIBUTION  OF  SNR  GAINS  FROM  AVERAGE  SUBARRAY  TO  ADJUSTED 
DELAY  ARRAY  BEAM  (1971  EVENTS;  WIDE -BAND  DATA) 


DELAY  ARRAY  BEAM  (1971  EVENTS; 


side  are  contributed  chiefly  by  events  v/hose  reference  sub¬ 
array  is  much  better  than  the  average  subarray.  Parenthe¬ 
tically,  it  was  demonstrated  in  Special  Report  No.  9  that  the 
degree  of  spread  in  subarray  beam  SNR  's  for  a  given  event 
varies  systematically  from  one  epicentral  region  to  another. 
Some  of  them,  however,  represent  events  with  relatively  poor 
signal  similarity  ~>etween  subarrays. 

Diversity-stack  beam  performance  relative  to  the  adjusted- 
delay  beam  is  1.  0  +  0.  1  dR  for  wide-band  data  and  slightly 
higher  for  filtered  data,  essentially  the  same  as  that  noted  in 
previous  reports.  Figure  11-5  shows  that  values  for  all  events 
are  distributed  closely  about  the  mean.  Negative  values  can 
be  explained  as  the  result  of  poor  coherence  among  the  stronger 
subarrays. 

Filter  improvement  is  consistently  more  variable  than  beam- 
fc  rming  improvement  as  a  direct  result  of  the  wide  variation 
ir  frequency  content  from  event  to  event  (Figures  II-6,  II-7). 

It  v\s  shown  in  Special  Report  No.  9  (197  3),  that  this  spectral 
variation  is  highly  dependent  on  source  region.  Standard  de¬ 
viations  are  on  the  order  of  4  dB  for  average  and  reference 
subarray  beam  and  both  array  beams.  Mean  gains  are  about 
8  dB  for  both  average  subarray  and  adjusted-delay  beams. 
Slightly  higher  gains  for  reference  subarray  and  diversity- 
stack  beams  are  due  to  decreased  attenuation  of  high  frequen¬ 
cies  in  the  signal. 

Mean  total  array  processing  gain,  (Figure  11-8)  is  18  +,4.  5  dB; 
including  the  7  dB  subarray  beamforming  gain  mentioned  be¬ 
fore  overall  SNR  gain  (from  unfiltered  singl  sensor  to  filtered 
array  beam)  is  about  25  ±  5  dB. 


DISTRIBUTION  OF  SNR  GAINS  FROM  ADJUSTED-DELAY  ARRAY  BEAM 
TO  DIVERSITY-STACK  BEAM  (1971  EVENTS;  WIDE-BAND  DATA) 


As  mentioned  previously,  all  the  above  numbers  refer  to  1971 
events,  for  which  "optimum"  processing  was  performed.  For  our  events 
from  1972,  we  obtained  an  adjusted-delay  array  beamforming  gain  of  7  _+ 4  dB 
relative  to  the  average  subarray  beam  for  filtered  signals.  Thus  we  estimate 
that  our  use  of  delays  from  nearby  master  events  results  in  a  decline  of  about 
3  dB  in  array  performance,  relative  to  the  ideal  case.  Obviously,  this  loss 
may  be  reduced  by  applying  interpolation  techniques  and  also  by  enlarging  the 
set  of  master  events.  Therefore,  our  statistics  from  1971  and  1972  may  con¬ 
veniently  be  viewed  as  upper  and  lower  bounds,  respectively,  of  operational 
array  performance. 

C.  PARTIAL  ARRAY  PERFORMANCE 

Performance  of  two  eight-subarray  partial  arrays  relative  to 
that  of  the  full  array  was  evaluated  on  the  basis  of  data  from  60  events.  Means 
and  standard  deviations  of  indicators  of  partial  array  performance  are  shown 
in  Table  II-2.  Figures  II- 9  and  II- 10  show  histograms  of  SNR  loss  from  full 
array  beam  to  partial  array  beam  for  the  two  partial  arrays.  Theoretical 
beamforming  improvement  by  a  partial  array  of  this  size  is  9.  0  dB,  a  de¬ 
crease  of  4.4  dB  from  the  13.4  dB  for  the  full  array. 

In  the  case  of  the  partial  array  composed  of  the  inside  rings 
of  the  full  array  (Figure  II- 1 ),  the  observed  decline  in  performance  averaged 
slightly  less  than  5  dB. 

The  second  partial  array  (Figure  II-2)  was  composed  of  sub¬ 
arrays  in  the  northeast  quadrant  of  the  array  which  have  shown  consistently 
higher  SNR's  than  the  rest  of  th  1  array  for  epicentral  regions  of  interest 
(Special  Report  No.  9,  1973).  The  average  performance  of  this  partial  array 
was  only  2  dB  worse  than  that  of  the  full  array. 

Mean  filter  gain  did  not  show  any  significant  differences  be¬ 
tween  the  two  partial  arrays  and  the  full  array. 


TABLE  II- 2 

NORSAR  SHORT  PERIOD  PARTIAL  ARRAY 
PERFORMANCE  STATISTICS 


Ratio  of  Partial  Array  Beam  SNR  to  Full  Array  Beam  SNR 

Rings  A  and  B 

Northeast  Quadrant 

Wide-Band 

-  4.  8  +  3.  2  dB 

-  2.  1  +  2.0  dB 

Standard  Filter 

-  4.  7  +  2.4  dB 

-  2.  0  +  1.  9  dB 

Adjusted-Delay  Beam  Filter  Gains 

Full  Array 

Rings  A  and  B 

Northeast  Quadrant 

8  .9  ±4.3  dB 

9.  1  +5,4  dB 

9.  1  ±  3.  6  dB 
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It  may  thus  be  concluded  that  a  partial  array  (of  8  subarrays) 
in  the  northeast  corner  of  the  NORSAR  array  will  have  a  detection  capability 
that  is  only  slightly  inferior  to  that  of  the  full  array  for  Eurasian  events. 

Since  the  aperture  of  the  partial  array  is  only  about  one-half  of  that  of  the  full 
array,  a  decline  in  the  event  location  accuracy  must  alio  be  expected.  Finally, 
it  should  be  st  essed  that  our  analysis  is  valid  for  Eurasian  events  only  In 
fact,  for  events  from  the  Western  Hemisphere,  the  northeast  quadrant  of  the 
array  no  longer  gives  the  highest  signal  amplitudes,  and  a  different  partial 
array  configuration  would  probably  be  optimum  for  that  case. 


11-18 


SECTION  III 

NORSAR  TELESEISMIC  DETECTION  CAPABILITY 


A.  INTRODUCTION 


One  of  the  main  objectives  of  this  study  was  to  determine  the 
detectability  of  P-waves  using  the  short-period  NORSAR  array.  For  this 
purpose  an  adjusted  delay  array  beam  was  formed  for  each  event,  and  the 
decision  detection/ no  detection  was  made  by  the  TI  analyst  after  visual  in¬ 
spection  of  the  signal  traces.  For  most  events  having  bodywave  magnitude 
lower  than  4.  5  prefiltering  of  the  signal  was  performed  using  our  standard 
filter,  which  is  similar  to  the  bandpass  filter  used  in  the  NORSAR  on-line 
Detection  Processor.  The  procedure  leading  to  the  selection  of  this  filter 
was  outlined  in  Special  Report  No.  6,  (1972).  The  NORSAR  SP  incremental 
detection  thresholds  were  then  estimated  by  observing  the  detection  percen¬ 
tages  as  a  function  of  event  magnitude. 


The  method  of  estimation  utilized  in  this  report  is  a  maximum 
likelihood  procedure  described  by  Ringdal  (1974).  Briefly,  this  method  as¬ 
sumes  that  the  probability  P(m)  of  a  seismic  station  detecting  an  event  of 
bodywave  magnitude  m  in  a  given  region  may  be  described  by  an  error 
function: 


m 


P(m)  =  (27t<t  ) 


j 


(t  -H) 
ZaZ 


dt 


(1) 


In  this  model,  then,  the  station  detection  potential  lor  events  from  a  specific 
region  is  characterized  by  the  corresponding  values  of  y,  and  cr  .  The  pro¬ 
blem  is  therefore  to  estimate  these  two  parameters. 


Ill- 1 


The  general  procedure  in  estimating  ^  and  a  is  as  follows: 

•  Obtain  a  reference  set  of  randomly  selected  events  of  various 
magnitudes. 

•  For  each  event  in  the  reference  set,  make  a  decision  as  to 
whether  or  not  the  station  has  detected  this  event. 

•  Establish  the  likelihood  function  for  the  observed  pattern  of 
decisions:  detection  versus  no  detection,  using  the  probability 
model  (1).  (This  is  easy  since  ail  decisions  can  be  assumed 
independent.  ) 

•  Find  the  values  of  the  parameters  that  maximize  the 

above  likelihood  function. 

When  applying  this  method  to  the  NORSAR  array,  some  con¬ 
sideration  must  first  be  given  to  the  reference  event  set.  It  is  essential  that 
the  events  in  the  reference  set,  for  any  given  magnitude,  represent  a  random¬ 
ly  chosen  subset  of  the  total  number  of  events  occurring.  In  this  way,  the 
percent  detected  will  be  an  unbiased  estimate  of  the  percentage  that  NORSAR 
would  detect  of  the  whole  event  population  for  each  magnitude. 

The  randomness  criterion  above  means  that  no  event  must  be 
chosen  with  any  a  priori  knowledge  as  to  whether  or  not  NORSAR  can  be  ex¬ 
pected  to  detect  this  event.  (Apart,  of  covrse,  from  magnitude  and  location 
information.  )  Thus  we  require  statistical  independence  between  NORSAR  and 
the  reporting  source  with  respect  to  event  detection  probability. 

In  order  to  achieve  this  independence,  we  deleted  all  events  in 
our  data  base  that  had  originally  been  selected  using  the  NORSAR  seismic 
bulletin  as  a  source.  We  also  excluded  events  reported  by  IS!/  where  NORSAR 
had  been  listed  as  the  primary  source,  (i.  e.  ,  events  that  would  probably  not 
have  been  reported  if  NORSAR  had  not  detected  them).  The  slight  bias  into- 
duced  by  the  fact  that  NORSAR  contributes  to  the  compilation  of  the  PDE 
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bulletins  was  considered  insignificant  for  our  purposes;  we  therefore  includ¬ 
ed  events  reported  by  PDE  in  our  reference  set, 

The  data  base  for  our  detection  threshold  estimation  was  in 
this  way  reduced  to  452  Eurasian  events.  17  3  events  were  located  on  the 
Kamchatka-Kuriles  arc;  38  of  these  were  not  detected  on  NORSAR  data.  Of 
the  remaining  2 79  events,  all  but  12  were  detected.  This  high  detection  rate 
is  due  partly  to  the  good  capability  of  the  NORSAR  array  for  Europe  and  Cen¬ 
tral  Asia,  partly  to  the  lack  of  low  magnitude  events  from  this  region  in  our 
data  base;  this,  again,  reflects  a  relatively  poor  coverage  provided  by  LASA. 
As  a  consequence,  it  has  not  been  feasible  to  perform  a  detection  threshold 
estimation  on  a  regional  basis  for  the  Eurasian  continent. 

For  large  events,  the  time  delay  adjustments  for  the  array 
beamforming  were  computed  by  cross-correlating  the  subarray  waveforms. 

For  all  low  magnitude  events,  time  delays  were  estimated  on  the  basis  of 
nearby  large  events.  Thus  it  may  be  assumed  that  it  would  have  been  possible 
to  obtain  more  precise  regional  corrections  if  a  larger  data  base  and  more 
elaborate  processing  methods  had  been  used.  As  a  consequence,  the  estimates 
of  the  NORSAR  short  period  detection  capabilities  presented  here  should  be  re¬ 
garded  as  conservative. 

Finally,  it  should  be  observed  that  all  our  estimates  are  in 
terms  of  LASA,  PDE,  or  ISM  magnitudes  v  lich  are  considered  to  be  mutu¬ 
ally  unbiased).  It  is  important  to  quote  the  source  stations  when  referring 
to  our  detection  threshold  estimates. 

B.  ESTIMATE  OF  THE  NORSAR  m,  DETECTION  THRESHOLD 

b 

Figure  III  - 1  presents  a  histogram  showing  the  NORSAR  detec¬ 
tion  performance  for  the  Kuriles-Kamchatka  arc  (Region  1)  together  with  the 
associated  maximum  likelihood  detection  curve.  The  90  percent  incremental 
threshold  estimate  is  around  mb=  4.  3.  The  epicentral  distances  from  NORSAR 
are  between  60  and  75  degrees  for  events  from  this  region. 
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FIGURE  III-  1 

NORSAR  SP  DETECTION  STATISTICS  FOR  THE 
KURILES-KAMCHATKA  REGION 
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A  similar  picture  for  Eurasia  apart  from  Kuriles-Kamchatka 
(Region  2)  is  given  in  Figure  III —2 .  Because  of  the  relatively  few  low  magni¬ 
tude  events  in  this  data  base,  the  detection  curve  cannot  be  estimated  with 
any  great  degree  of  confidence  for  low  values.  However,  a  reasonably 
accurate  estimate  of  the  90  percent  incremental  detection  threshold  may  still 

be  obtained;  it  appears  to  lie  around  m  =  4,  0.  Most  of  our  events  from 

b 

Region  2  are  from  Europe  and  Central  Asia,  with  epicentral  distances  to 
NORSAR  ranging  from  20  to  50  degrees. 

It  shoulc  be  pointed  out  that  several  of  our  events  for  Region  2 

have  relatively  unreliable  mfa  estimates.  This  is  due  to  the  fact  that  PDE 

includes  near  regional  stations  in  their  mb  computation,  and  therefore  will 

often  report  too  high  a  magnitude.  In  particular  this  applies  to  events  for 

which  only  a  few  stations  report  an  amplitude,  since  those  stations  are  then 

likely  to  have  a  favorable  radiation  pattern  for  the  event  in  question.  As  a 

partial  solution  to  this  problem,  we  decided  to  delete  events  with  m  values 

b 

based  on  just  one  PDE  station  from  our  maximum  likelihood  estimation.  Thus 
one  event  of  mb  =  4.  5  and  one  of  mb  =  4.  8  (IT A/ 035/ 04N)  that  were  not  de¬ 
tected  were  eliminated  from  consideration. 

The  detection  statistics  for  all  our  events  from  Eurasia  arc 
presented  in  Figure  111-3.  The  90  percent  incremental  detection  threshold  is 
close  to  mb  =  4.2,  and  can  be  given  with  good  confidence. 

In  order  to  examine  a  possible  seasonal  influence  on  the  NORSAR 
detection  capability,  the  data  base  was  split  into  a  summer  event  suite  (origin 
date  April  through  September)  and  a  winter  event  suite  (origin  date  October 
through  March).  The  corresponding  detection  curves  are  presented  in  Figures 
II1-4  and  111-5.  The  incremental  90  percent  threshold  estimates  are  between 
4.  1  and  4.  2  for  the  summer  event  suite  and  between  4.  2  and  4.  3  for  the  winter 
events.  This  difference  of  about  0.  1  mb  units  may  be  attributed  to  the  gen¬ 
erally  lower  seismic  noise  level  during  the  summer  months. 
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It  is  interesting  to  notice  that  the  difference  found  here  is  less 

than  the  0.  3  m,  units  of  seasonal  difference  observed  for  the  NORSAR  LP 
b 

detection  thresholds  (Special  Report  No.  7,  1973).  This  confirms  the  obser¬ 
vation  from  our  SP  noise  analysis  (Special  Report  No.  6,  1972)  that  the  seis¬ 
mic  noise  filtered  by  the  standard  filter  exhibits  less  seasonal  variation  than 
the  wide-band  noise  level.  This  again  is  due  to  a  shift  in  the  microseismic 
peak  towards  lower  frequencies  as  the  noise  level  increases. 

C.  ESTIMATE  OF  THE  NORSAR  OPERATIONAL  DETECTION 
THRESHOLD 

A  comparison  between  the  detection  results  from  Tl's  analysis 
of  NORSAR  data  and  the  NORSAR  seismic  bulletin  was  presented  in  Special 
Report  No.  9i  (1973).  For  the  time  period  covered  by  that  investigation 
(January  to  March  1972)  it  appeared  that  NORSAR  was  operating  well  below 
its  potential.  Considering  that  the  weekly  number  of  events  reported  by 
NORSAR  has  increased  significantly  since  then,  it  would  be  of  interest  to  de¬ 
termine  the  present  operational  detection  level  for  the  array. 

For  this  purpose,  we  analyzed  the  NORSAR  performance  for 

two  event  swarms,  one  from  South  Honshu,  December  3-20,  1972;  the  other 

from  the  Kurile  Islands,  June  17-30,  1973.  LASA  detections  and  m,  values 

b 

were  used  as  a  reference,  and  detections  occurring  when  NORSAR  was  out  of 
operation  were  deleted.  In  this  way  statistics  were  compiled  concerning 
detection/no  detection  information  as  presented  in  Table  III  - 1  and  Figures 
III- 6  and  II1-7. 

One  interesting  observation  from  Table  111-1  is  that  NORSAR 
detected  more  events  than  LASA  from  the  Kuriles  swarm,  while  this  picture 
was  reversed  for  the  South  Honshu  aftershocks.  Both  of  these  swarms  were 
at  almost  identical  epicentral  distance  from  LASA  and  NORSAR.  The  main 
explanation  is  the  NORSAR  noise  level,  which  was  slightly  less  than  0.  2  m/i 
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TABLE  111-  1 


NORSAR  AND  LASA  EVENT  DETECTION  PERFORMANCE 
FOR  EARTHQUAKE  SWARMS  FROM  SOUTH  HONSHU 
(DECEMBER  3-20,  1972)  AND  THE  KURILE 
ISLANDS  (JUNE  17-30,  1973) 


Honshu 

Swa  rm 

Kurile  s 

Swa  rm 

Distance  from  NORSAR  (deg) 

78 

70 

Distance  from  LASA  (deg) 

80 

69 

LASA  Total  Detected  Events 

192 

364 

NORSAR  Total  Detected  Events 

133 

452 

Common  Events 

106 

284 
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RMS  on  the  bandpass  filtered  (1.2  -  3.2  Hz)  array  beam  in  December  1972, 
while  less  than  0.  1  m /i  in  the  last  part  of  June  1973.  This,  of  course, 
caused  the  90  percent  detection  level  for  the  Kuriles  swarm  to  be  significantly 
lower  than  that  of  the  Honshu  swarm,  (4.  1  compared  to  4.  3).  The  difference 
is  even  greater  in  the  50  percent  thresholds  (3.  5  and  3.9). 

It  is  interesting  to  notice  the  much  larger  spread  in  the  detec¬ 
tion  curve  for  the  Kuriles  swarm  compared  to  the  Honshu  curve,  (rr=  0.47  to 
(T~  0.  31).  Part  of  this  difference  may  be  due  to  a  greater  variation  in  the 
seismic  noise  level  during  June  1973,  but  mostly  it  seems  to  reflect  more 
variability  in  the  source  mechanisms  for  the  Kuriles  aftershocks. 

In  order  to  compare  Tl's  detectability  estimates  to  the  observed 
array  performance,  it  is  most  correct  to  compare  the  50  percent  detectability 
limits.  This  is  because  there  is  an  inherently  larger  spread  in  the  detect¬ 
ability  curve  when  the  events  span  z  larger  region  and  a  longer  period  of  time 
as  opposed  to  an  aftershock  sequence.  By  compensating  for  the  differences 
between  the  actual  noise  levels  during  the  two  event  swarms  and  the  annual 
average  level  of  0.  16  m^  on  the  filtered  array  beam  (Special  Report  No.  6, 
1972)  we  obtained  estimates  of  3.7  and  3.  85  for  the  50  percent  operational 
detection  levels  for  Kuriles  and  Honshu,  respectively,  during  normal  noise 
conditions.  These  numbers  are  only  slightly  higher  than  the  TI  estimate  of 
3.  65  for  the  Kuriles-Kamchatka  arc.  Thus  it  appears  that  NORSAR  presently 
operates  very  close  to  its  optimum  capability  for  this  region. 

It  has  not  been  possible  to  determine  the  present  NORSAR 
operational  capability  for  near  regional  events  in  a  way  similar  to  the  above 
procedure.  Our  investigations  presented  in  Special  Report  No.  9  (1973) 
showed  that  the  NORSAR  performance  was  significantly  below  the  array's 
potential  for  close-in  regions  at  that  time.  However,  it  is  believed  that  im¬ 
proved  regional  time  delay  corrections  and  the  implementation  of  envelope 
beamforming  in  the  NORSAR  Detection  Processor  may  have  altered  the  pic¬ 
ture  since  then. 


One  final  remark  seems  highly  relevant  in  view  of  the  observed 
array  detection  capabilities  for  the  two  earthquake  swarms  examined  in  this 
section:  The  v.  riability  of  the  NORSAR  seismic  noise  level  is  considerable, 
even  within  the  shori  period  detection  filter  band,  and  no  general  s*ateme'  t 
about  the  NORSAR  detection  capability  should  be  made  without  taking  due  re¬ 
servations  for  this  factor. 


SECTION  IV 

SHORT  PERIOD  DISCRIMINATION 

A.  INTRODUCTION 

This  section  presents  the  results  achieved  by  applying  five 
standard  short  period  discriminants  to  NORSAR  data.  The  discriminants  are 
briefly  described  as  follows: 

1.  P30  Mean  Square 

This  discriminant,  which  is  a  measure  of  event  complexity,  is 
computed  by  crosscorrelating  4  seconds  of  the  waveform  (beginning  a  few  points 
before  P-wave  onset)  with  the  next  30  seconds  of  the  waveform  and  with  the 
noise  preceding  the  signal.  A  mean  square,  weighted  by  the  lag,  is  then  com¬ 
puted  from  the  correlations  over  both  30  seconds  of  the  noise  and  30  seconds 
of  the  signal.  The  noise  mean  square  is  subtracted  from  the  signal  mean 
square  to  obtain  the  discriminant  used  (Texas  Instruments  Incorporated,  1971). 

2.  Autocorrelation  Mean  Square 

This  discriminant  is  also  a  measure  oi  complexity.  The  auto¬ 
correlations  of  a  30-second  noise  gate  and  of  a  30-second  signal  gate  are  com¬ 
puted  and  a  weighted  mean  square  then  derived  from  these  correlations  for  the 
noise  and  signal.  The  discriminant  is  derived  from  the  signal  mean  square 
minus  the  noise  mean. 

3.  Envelope  Difference 

This  discriminant  is  also  derived  from  the  P30  correlation  by 
computing  the  mean-square  difference  between  the  envelope  correlation  and  a 
fixed  decaying  exponential,  the  decay  rate  of  which  is  the  average  rate  for  an 
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ensemble  of  16  explosions  recorded  at  LASA.  As  with  the  first  two  statistics, 
envelope  difference  is  a  measure  of  complexity. 

4.  Dominant  Period 

This  discriminant  is  computed  by  finding  the  cycle  in  the  wave¬ 
form  with  a  maximum  absolute  amplitude;  the  dominant  period  is  the  duration 
of  this  cycle  in  seconds.  This  parameter  can  be  estimated  with  come  confidence, 
even  for  events  with  a  relatively  low  signal-to-noise  ratio.  The  dominant  per¬ 
iod  discriminant  is  a  rough  measure  of  Bpectral  energy  distribution. 

5.  Spectral  Ratio 

This  discriminant  is  derived  from  the  signal  power  soectrum 
over  a  gate  beginning  just  before  the  signal  arrival.  The  power  spectrum  is 
smoothed  over  three  frequency  points,  and  the  power  in  three  bands  is  com¬ 
puted;  Band  1  ;  0.  -0.  55  Hz;  Band  2  :  0.  55  -  1. 5  Hz;  Band  3  ;  1.  5  -  5.  0  Hz. 

These  bands  have  been  selected  based  on  NORSAR  data.  Spectral  ratios  com¬ 
puted  were  Band  3  to  Band  2  and  Band  3  to  Band  1,  respectively. 

In  order  to  evaluate  the  performance  of  the  individual  discri¬ 
minants  as  well  as  the  possible  combined  criteria,  a  measure  of  the  separa¬ 
tion  achieved  between  earthquakes  and  presumed  explosions  was  computed 
for  each  discriminant.  Since  all  our  discriminants  are  two-dimensional, 
(discriminant  value  versus  ) ,  the  problem  of  obtaining  such  a  measure 
reduces  to  measuring  the  separation  of  two  point  sets  in  a  plane.  The  follow¬ 
ing  method  was  adapted  (see  Figure  IV-1). 

1.  For  any  given  straight  line,  the  distances  from  the  line  were 
computed  for  all  points  corresponding  to  events  in  the  earthquake  and  pre¬ 
sumed  explosion  populations, 

2.  The  two  sets  of  real  numbers  thus  obtained  were  considered  as 
sampled  from  two  Gaussian  populations;  one  N^,  (earthquakes)  and  one 
N(M2'  (Presumed  explosions). 
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FIGURE  IV -1 

ILLUSTRATION  OF  A  METHOD  TO  MEASURE  THE  SEPARATION 
BETWEEN  TWO  POINT  SETS  IN  A  PLANE 


3.  The  straight  line  that  gave  two  Gaussian  populations  with  the 

best  separation  was  chosen,  and  the  measure  of  separation  was  defined  as 
the  corresponding  probability  of  correct  classification. 


There  is  no  unique  method  of  computing  the  separation  of  two 
Gaussian  distributions  with  unequal  variances.  We  chose  under  3  above  to 
measure  the  separation  by  computing  the  worst  case  probability  of  misclass- 
ification,  applying  the  well-known  minimax  principle  (see  e.  g.  van  Trees 
(1968)).  The  minimax  separation  threshold  between  the  two  populations  is 


given  as 


tMlCT2 
<7,  + 


and  the  associated  probability  of  error  is 

P(miss)  =  P  (false  alarm)  =  <I> 

where  <J>  is  the  standard  cumulative  Gaussian  distribution  function. 

Thus  it  is  seen  that  for  the  above  threshold  T  ,  the  conditional 
probabilities  of  correctly  identifying  an  earthquake  and  correctly  identifying 
a  presumed  explosion  are  identical,  and  therefore  provide  a  well-defined 
measure  of  the  separation  of  the  two  populations. 

Finally,  it  may  be  observed  that  the  above  method  suggests  an 
easy  way  to  evaluate  the  separation  properites  of  multivariate  discriminants. 
In  fact,  since  the  method  reduces  a  two-dimensional  discriminant  to  two  sets 
of  real  numbers,  it  can  be  applied  successively  to  reduce  any  given  array  of 
discriminant  values  to  the  one-dimensional  case. 


B.  NORSAR  SHORT  PERIOD  DISCRIMINATION  RESULTS 

Short  period  discriminant  values  for  the  discriminants  defined 
in  Subsection  IV- A  are  plotted  as  a  function  of  body-wave  magnitude  for  a 
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total  of  414  events  in  Figures  IV-2  through  IV- 13.  Shallow  earthquakes  and 
earthquakes  of  unknown  depth  are  represented  by  a  cross.  Deep  earthquakes 
(of  depths  greater  than  100  km)  are  denoted  by  a  triangle.  Presumed  explo¬ 
sions  are  indicated  by  an  asterisk.  Events  from  the  Western  Hemisphere 
are  surrounded  by  a  circle. 

The  "best  separation  line"  as  defined  in  Subsection  1V-A  be¬ 
tween  shallow  Eurasian  earthquakes  and  Eurasian  presumed  explosions  is 
drawn  for  each  of  the  Figures  IV-2  through  IV- 13.  Note  that  this  line  has 

been  found  on  the  basis  of  all  events  with  m  ^  4.4;  this  is  because  all  our 

b 

presumed  explosions  are  of  at  i^ast  this  magnitude. 

Table  IV- 1  lists  the  results  obtained  when  evaluating  the  var¬ 
ious  criteria  according  to  the  method  introduced  in  Subsection  1V-A.  The 
"obabilities  of  correct  identification  are  generally  between  80  and  90  percent. 
For  all  discriminants  except  the  envelope  difference,  better  separation  is 
achieved  on  the  array  beam  than  when  using  the  reference  subarray,  although 
the  differences  are  not  large.  As  can  be  expected,  the  values  of  the  discri¬ 
minants  at  the  subarray  and  array  level  are  highly  correlated  with  correlation 
coefficients  up  to  0.  90.  The  spectral  ratio  criteria  show  the  best  performance 
both  on  the  subarray  and  array  level,  with  Band  3/ Band  1  appearing  to  give  the 
best  separation.  This  particular  criterion  was  found  to  identify  correctly  92 
percent  of  the  events.  However,  some  reservations  should  be  made  on  the 
grounds  that  this  spectral  ratio  to  a  large  degree  measures  signal-to-noise 
ratio.  The  good  performance  may  thus  in  part  be  due  to  the  short  epicentral 
distance  to  our  presumed  explosions,  which  again  causes  a  high  SNR  relative 
to  m^  for  these  events  (Figure  IV-14). 

Table  IV-2  presents  the  separation  properties  achieved  by  com¬ 
bining  any  two  of  our  criteria  at  the  array  beam  level.  The  improvements  are 
not  large;  the  best  performance  (93  percent  correct  identification)  is  obtained 
by  combining  Log  P30  with  Spectral  Ratio  3/1.  Incidentally,  the  separation 
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REFERENCE  SUBARRAY  P30  MEAN  SQUARE  DISCRIMINANT 
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FIGURE  IV -3 

ADJUSTED-DEDAY  BEAM  P30  MEAN  SQUARE  DISCRIMINANT 
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FIGURE  IV -4 

REFERENCE  AUTOCORRELATION  MEAN  SQUARE 
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FIGURE  IV -5 

ADJUSTED- DELAY  BEAM  AUTOCORRELATION  MEAN  SQUARE 
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FIGURE  IV -6 

REFERENCE  ENVELOPE  DIFFERENCE  (COUNTS  SQUARED) 
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FIGURE  IV -8 

REFERENCE  DOMINANT  PERIOD  (SECONDS) 
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FIGURE  IV -9 

ADJUSTED-DELAY  BEAM  DOMINANT  PERIOD  (SECONDS) 
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FIGURE  IV-10 

REFERENCE  SPECTRAL  RATIO  [(1.5  -  5)/(0  -  0.  55)  J  (dB) 
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FIGURE  IV -12 

REFERENCE  SPECTRAL  RATIO  [(1.5-  5)/(0.  55  -1.5)]  (dB) 
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TABLE  IV-1 

PROBABILITIES  OF  CORRECT  IDENTIFICATION  (SHALLOW 

EARTHQUAKE  vs.  PRESUMED  EXPLOSION  USING 

SHORT  PERIOD  DISCRIMINANTS  FOR 

EURASIAN  EVENTS  OF  m.  >  4.  4 

b 


Reference 

Subarray 

Adj.  Delay 
Array  Beam 

Correlation 
Between  Ref.  SA 
And  AB  Criteria 

Log  P30 

0.  831 

0.  861 

0.  86 

Log  Autocorrelation 

0.  816 

0.  847 

0.  66 

Envelope  Difference 

0.  796 

0.  784 

0.  90 

Dominant  Period 

0.  816 

0.  878 

0.  66 

Spectral  Ratio  3/1 

0.  874 

0.  917 

0.  84 

Spectral  Ratio  3/2 

0.  875 

0.  898 

0.  89 

All  Discriminants  Combined 

'  “ —  -  .1 

0.  897 

0.  934 
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FIGURE  IV - 14 

SIGNAL-  TO -NOISE  RATIO  VERSUS  MAGNITUDE  (m  )  FOR 
WIDE-BAND  ARRAY  BEAMS  AT  NORSAR  b 
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TABLE  IV- 2 

PROBABILITIES  OF  CORRECT  IDENTIFICATION  (SHALLOW  EARTHQUAKE 
vs.  PRESUMED  EXPLOSION)  USING  COMBINATIONS  OF  TWO 
SHORT  PERIOD  DISCRIMINANTS  APPLIED 
AT  THE  ARRAY  BEAM  LEVEL 
(EURASIAN  EVENTS  OF  m  >4.  4) 


Combined  With: 

I 

2 

3 

4 

5 

6 

1. 

Lop  P30 

0.  861 

0.  882 

0.  862 

0.  910 

0.  932 

0.  929 

2. 

Lop  Autocorrelation 

0.  847 

0.  857 

0.  892 

0.  922 

0.  908 

3. 

Envelope  Difference 

0.  857 

0.  784 

0.  885 

0.  923 

0.  906 

4. 

Dominant  Period 

0.  885 

0.  878 

0.  922 

0.  906 

5. 

Spectral  Ratio  3/1 

0.  922 

0.  923 

0.  922 

0.  917 

0.  927 

6. 

Spectral  Ratio  3/2 

0.  908 

0.  906 

0.  906 

0.  927 

0.  898 

I  V-  20 


provided  by  this  combination  is  essentially  equal  to  that  obtained  by  combining 
all  the  criteria  (Table  IV  —  1 ) ;  this  shows  that  there  is  not  much  more  to  be 
gained  by  any  multivariate  discriminant  with  our  definition  of  separation. 

In  order  to  measure  the  interdependence  of  our  SP  discrimin¬ 
ants,  we  computed  the  correlation  coefficients  between  each  pair  of  discri¬ 
minants  at  the  array  beam  level  as  shown  in  Table  IV-3.  Presumed  explosions 
and  earthquakes  were  treated  separately.  Not  unexpectedly  these  computations 
show  that  the  three  discriminants  based  on  spectral  characteristic  s  are  strong¬ 
ly  interrelated  (correlation  coefficients  0.  6  to  0.  8),  as  is  also  the  case,  to  a 
somewhat  lesser  extent,  with  the  complexity  criteria.  However,  correlations 
across  these  two  classes  arc  generally  low,  although  in  most  cases  positive. 
This  confirms  the  observation  from  Table  IV-2  that  the  highest  improvements 
in  separations  are  generally  obtained  by  combining  one  spectral  and  one  com¬ 
plexity  type  discriminant. 

As  an  alternative  way  of  combining  complexity  with  spectral 
information,  it  was  attempted  to  apply  the  standard  filter  prior  to  computation 
of  the  complexity  discriminants.  This  resulted  in  some  improvement  in  sep¬ 
aration,  but  the  performances  were  still  not  comparable  to  those  of  the  spectral 
c  riter  ia. 

Finally,  several  frequency  bands  were  examined  to  see  if  alter¬ 
native  choices  could  improve  appreciably  the  spectral  ratio  discriminants.  It 
was  found,  however,  that  it  was  not  possible  to  obtain  any  significant  improve¬ 
ment  this  way,  and  we  therefore  may  consider  the  measured  performance  of 
our  spectral  ratios  as  typical  for  this  kind  of  discriminant. 

As  was  stated  in  Section  I,  the  combined  performance  of  short 
and  long  period  discriminants  at  NORSAR  will  be  evaluated  in  a  forthcoming 
report.  'I  hat  report  will  also  provide  detailed  case  studies  of  events  that  fail 
to  tie  classified  properly  by  the  combined  criteria. 
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TABLE  IV- 3 


CORRELATION  COEFFICIENTS  BETWEEN  SP  ARRAY  BEAM 
DISCRIMINANT  VALUES  FOR  EURASIAN  EVENTS  OE  mi  >4.4. 
NUMBERS  BELOW  THE  DIAGONAL  REPRESENT  THE  SHALLOW 
EARTHQUAKE  POPULATION.  WHILE  NUMBERS  ABOVE  THE 
DIAGONAL  CORRESPOND  TO  PRESUMED  EXPLOSIONS. 


Correlated  With 


I 

2 

3 

4 

5 

6 

L 

Log  P30 

1.  00 

0.  59 

0.  80 

0.  26 

0.  31 

0.  15 

2. 

Log  Autocorrelation 

0.  27 

1.  00 

0.  49 

0.  43 

0.  49 

0.  37 

3. 

Envelope  Difference 

0.  71 

0.  14 

1.  00 

0.  23 

-0.  03 

0.  18 

4. 

Dominant  Period 

0.  20 

0.  47 

0.  17 

1.  00 

0.  63 

0.  85 

5. 

Spectral  Ratio  3/1 

0.  25 

0.  40 

0.  15 

0.  68 

1.  00 

0.  56 

6. 

Spectral  Ratio  3/2 

0.  04 

0.  45 

0.  08 

0.  61 

0.  74 

I.  00 

SECTION  V 

SUMMARY  AND  CONCLUSIONS 


This  section  summarizes  the  results  achieved  during  the  two 
and  one-half  year  NORSAR  short  period  evaluation  program. 

1 .  Data  Quality 

The  SP  data  recorded  at  NORSAR  has  been  of  consistently  high 
quality  throughout  the  evaluation  period,  which  spanned  the  time  interval 
from  March  1971  through  December  1972.  On  the  basis  of  more  than  500 
events  and  70  noise  samples  processed  by  TI  for  this  period,  the  following 
observations  were  made: 

•  Data  was  available  from  NORSAR  for  more  than  95  percent  of 
the  tim-  intervals  requested  by  TI. 

•  In  most  cases  at  least  20  subarrays  were  operational.  The 
worst  data  loss  for  a  single  event  was  33  out  of  132  sensors. 

•  Data  spikes  were  observed  for  10  events,  but  these  events 
could  still  be  processed. 

•  Phase  reversals  affected  8  sensors  during  parts  of  1971,  but 
was  not  observed  on  1972  data. 

•  The  SP  seismometers  appeared  to  be  well  equalized  across  the 
NORSAR  array. 

2.  Noise  Analysis 

The  following  conclusions  concerning  noise  analysis  are  based 
on  72  sample  intervals: 
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•  The  noise  spectral  shape  is  very  simple,  with  a  peak  at  about 

3  to  6  seconds  and  a  rapid  fall-off  toward  shorter  periods.  The 
spectral  shape  does  not  change  significantly  across  the  array. 

•  Noise  levels  are  very  similar  across  the  array.  Maximum 
single  sensor  variations  typically  are  +_  6  dB,  and  most  sensors 
are  within  3  dB  of  the  average  single  sensor  level.  Variation 
among  subarray  beam  noise  levels  is  4^2  dB. 

•  Wideband  RMS  noise  level  shows  a  significant  variation  with 
time,  and  correlates  strongly  with  storm  activity  in  the  North 
Atlantic  Ocean.  The  spectral  peak  generally  shifts  towards 
lower  frequencies  as  the  noise  level  increases,  Wintertime 
wideband  noise  levels  are  on  the  average  6  dB  higher  than 
summertime  levels;  this  difference  is  less  evident  when  the 
"standard"  bandpass  filter  is  applied. 

•  Typical  RMS  noise  levels  are:  0.  5  mfjL  ±  6  dB  for  the  wide  band 
array  beam.  0.  12  m/i  4.  3  dB  for  the  array  beam  through  the 
standard  filter.  This  last  number  is  about  a  factor  of  2  higher 
than  the  detection  band  noise  level  for  LASA. 

•  Multiple  coherence  levels  within  a  subarray  are  low  except  at 
the  3  to  6  seconds  microscismic  peak.  Inter-subarray  multiple 
coherencies  are  low  over  the  entire  0  to  5.  0  Hz  band. 

3.  Signal  Analysis 

Our  conclusions  from  the  signal  analysis  are  based  upon  the 
processing  of  S67  events;  and  can  be  summarized  as  follows: 

4'  Except  for  a  few  close-in,  high  frequency  events  signal  simil¬ 

arity  is  good  within  a  subarray.  Among  subarrays,  however, 
similarity  is  quite  variable. 
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Amplitude  variations  across  the  array  are  large,  typically  4:1, 
whik  variations  as  high  as  10:1  have  been  observed  for  Kazakh 
events.  The  amplitude  patterns  are  strongly  dependent  upon 
source  ocation,  but  consistent  behavior  is  generally  seen  with¬ 
in  narrow  regions.  It  appears  that  most  of  the  amplitude  vari¬ 
ations  may  be  explained  by  scattering  effects  due  to  the  irregular 
structure  of  the  Mohorovicic  discontinuity  underneath  the  NORSAR 
array. 

Tinn  delay  anomalies  (deviation  from  plane  wave  propagation 
along  the  great  circle  path)  are  not  significant  for  subarray 
beamforming.  Anomalies  are  significant,  however,  between 
subarrays  and  are  occasionally  as  large  as  1  second.  Consis¬ 
tent  s.-ts  of  anomalies  can  in  general  be  obtained  for  all  regions 
except  those  within  30°  cpicentral  distance  of  NORSAR. 

Fime -domain  signal  traces  from  various  regions  show  as  ex¬ 
pected,  a  general  tendency  towards  lower  complexity  as  the 
epicentr al  distance  (A)  increases.  Exceptions  to  this  rule 
were  some  high  complexity  signals  observed  for  Kamchatka 
events  (  \=  65°)  and  Taiwan  events  (  \=  80°). 

Signal  spectral  characteristics  show  strong  regional  variations, 
even  between  regions  very  close  together,  and  do  not  always 
follow  the  expected  tendency  towards  lower  frequencies  as  the 
epicentr  al  distance  increases.  Significant  high  frequency 
energy  (2  Hz  or  more)  was  observed  for  events  from  Greece 
(  \  =  25  ),  1  adzhik  (  \  -  40  )  and  the  Kurile  Islands  (  \  =  70°). 
Signals  of  dominant  low  frequency  (lower  than  1  Hz)  were  seen 
mainly  from  Italy  (  \=  20°),  Turkey  (  \=  25°),  Kirgiz  (  \=  45°) 
and  Taiwan  (  \  =  80°). 
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Our  limited  ensemble  of  Western  Hemisphere  events  show 
substantially  less  high  frequency  energy  than  the  Eurasian 
events. 


NORSAR  oody  wave  magnitudes  average  about  0.  2  units 
lower  than  either  PDE  or  LASA  values,  with  a  standard  devia¬ 
tion  of  0.  3  around  this  bias.  It  appears  that  this  negative  bias 
may  be  explained  as  signal  loss  in  array  beamforming.  The 
PDE-NORSAR  differences  appear  to  be  larger  at  low  mag¬ 
nitudes;  this  is  believed  to  be  because  PDE  in  those  cases  com¬ 
putes  an  m  based  upon  only  a  few  stations  with  favorable  ra¬ 
diation  patterns,  thus  resulting  in  too  high  m  values. 

b 

Array  Processing  Performance 


vN  noise  rejection  is  achieved  over  the  entire  0  to  5  Hz  band 
both  for  subarray  and  array  beamforming.  The  only  exceptions 
to  this  rule  are  occasional  strong  Rayleigh  wave  noise  fields 
(3-6  second  periods)  during  storm  activity  in  the  North  Atlantic 
Ocean  that  show  coherency  at  the  subarray  level.  Thus  noise 
rejection  totals  about  21  dB  (8  for  subarray  and  13  for  array 
beamforming). 


Signal  d 'gradation  for  subarray  beamforming  is  1  dB  +  0.  5  dB 
in  the  detection  frequency  band. 

Signal  degradation  for  array  beamforming  is  quite  variable, 
but  in  the  teleseismic  zone  the  following  values  were  found: 

3  dB  4;  2  dB  for  wide  band  signals 

3.  S  dB  ±2  dB  in  the  detection  frequency  band. 

Diversity-stack  beamforming  gives  the  following  SNR  im¬ 
provement  over  the  adjusted -delay  array  beam: 
i.  0  dB  ±0.9  dB  for  wide  band  signals 
1.6  dB  ±  1 . 0  dB  in  the  detection  frequency  band. 
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•  For  detection  of  Eurasian  events,  a  filter  with  corner  fre¬ 
quencies  at  about  1.2  and  2.  8  Hz  and  a  very  sharp  rolloff  at 
low  frequencies  appears  to  be  about  optimum.  This  "standard" 
filter  is  similar  to  the  1.2-3.  2  Hz  bandpass  filter  used  in  the 
NOKSAR  on-line  Detection  Processor. 

•  Gain  in  SNR  from  applying  the  standard  filter  was  highly  vari¬ 
able  and  showed  as  expected  a  strong  regional  dependence. 
Average  value  was  found  to  be  8+4  dB  both  at  the  subarray 
and  array  beam  levels. 

•  The  total  net  gain  of  the  NORSAR  array;  i.  e.  ,  the  SNR  im¬ 
provement  from  the  average  wide -band  single  sensor  to  the 
adjusted-delay  array  beam  filtered  with  the  standard  filter 
was  found  to  be  25  dB  +.  5  dB. 

•  The  performance  of  two  partial  NORSAR  arrays,  each  sonsist- 
ing  of  eight  subarrays,  were  evaluated  by  examining  SNR  losses 
for  60  Eurasian  events  relative  to  the  full  array.  A  partial 
array  consisting  of  the  A  and  B  rings  gave  an  average  loss  of 

4.  7  dB  SNR,  while  a  partial  array  situated  in  the  Northeast 
corner  of  NORSAR  averaged  a  loss  of  only  2  dB. 

5.  Event  Detection  Capabilities 

Event  detection  thresholds  were  estimated  on  the  basis  of  452 
processed  events  from  1971  and  1972  that  had  been  reported  by  sources  in¬ 
dependent  of  NORSAR.  A  maximum-likelihood  method  was  utilized  in  the 
estimation  procedure.  The  following  results  were  obtained;  (Note  that  all 
the  threshold  estimates  arc  relative  to  DDE,  LASA,  or  ISM  magnitudes.  ) 


V  -5 


•  90  percent  incremental  mfa  detection  threshold  for  all  of 
Eurasia  combined  is  approximately  4.2.  This  conforms  well 
to  the  corresponding  level  of  3.  9  for  LASA  reported  by  Dean 
et  al.  ,  (1971),  considering  that  the  NORSAR  noise  level  in  the 
detection  band  is  about  a  factor  of  2  higher  (0.  3  units). 

•  For  the  Kuriles -Kamchatka  arc  (epicentral  distance  60-70 
degrees)  the  90  percent  threshold  is  slightly  below  4.  3.  The 
average  value  for  the  remainder  of  Eurasia  (distances  generally 
20-59  degrees)  is  around  4.  0. 

•  The  winter  90  percent  threshold  was  found  to  be  slightly  higher 
(0.  1  m^  units)  than  the  summer  level.  This  difference  is  at¬ 
tributed  to  seasonal  variations  in  the  seismic  noise  level. 

•  A  theoretical  estimate  of  the  NORSAR  detection  threshold 
based  upon  seismic  noise  levels  and  measured  processing 
Josses  gave  results  consistent  with  the  direct  method. 

•  The  operational  event  reporting  performance  of  the  NORSAR 
system  was  found  to  be  well  below  array  capability  in  early 
1972,  especially  for  near  regional  events.  A  significant  im¬ 
provement  was  observed  one  year  later  for  the  Japan-Kuriles 
region,  where  the  array  then  appeared  to  be  operating  at  close 
to  optimum  capacity. 

6.  Short  Period  Disc  r  imination 

Five  standard  short  period  discriminants  were  applied  to  a 
total  of  414  events,  including  31  presumed  explosions,  27  of  which  were  from 
Eurasia.  The  main  results  are  as  follows: 

•  Our  SP  discriminants  do  not  appear  to  work  well  for  events 
from  the  Western  Hemisphere. 
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Discriminants  based  on  spectral  energy  distribution  seem  to 
be  superior  to  discriminants  based  upon  the  complexity  of  the 
signal  waveform. 

No  single  discriminant  was  able  to  separate  completely  be¬ 
tween  presumed  explosions  and  earthquakes.  The  best  separa¬ 
tion  was  obtained  by  considering  the  spectral  ratio  of  energy 
in  the  bands  1.5  -  5.  0  Hz  and  0-0.  55  Hz,  although  reserva¬ 
tions  must  be  taken  due  to  a  possible  bias  caused  by  the  high 
SNR  values  for  all  events  in  the  presumed  explosion  population. 

A  combination  of  SP  criteria  yielded  some  improvement  in 
separation,  but  no  substantial  change.  The  best  improvement 
was  obtained  by  combining  one  complexity  discriminant  with 
one  spectral  discriminant. 

A  preliminary  study  of  the  performance  of  short  period  dis¬ 
criminants  versus  that  of  M  -  m,  and  other  SP-LP  criteria 

s  b 

gave  the  expected  result  that  the  latter  ones  in  general  pro¬ 
duce  a  better  separation  between  earthquakes  and  presumed 
explos  ions. 
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APPENDIX  A 


THE  COMBINED  1971  -  1972  DATA  BASE 


All  events  processed  during  the  NORSAR  short-period  evalua¬ 
tion  are  listed  in  Table  A-l.  The  events  are  in  chronological  order  starting 
on  March  1,  1971  and  ending  on  December  28,  1972.  The  parameters  listed 
are:  event  designation,  date,  origin  time  (GMT),  latitude  and  longitude, 
depth  (km)  and  body-wave  magnitude  as  reported  by  the  source  institution. 

The  following  abbreviations  are  used  for  the  original  source  of 

info  rmation: 

P  :  Parameters  taken  from  Preliminary  Dete rmination  of  Epicenters 

Monthly  Summary 

I  :  Parameters  taken  from  the  International  Seismic  Month  event 

list 

S  :  Parameters  taken  from  the  SDAC/ L,ASA  Weekly  Summary 

N  :  Parameters  taken  from  the  NORSAR  Seismic  Event  Summary 

Note  that  in  the  cases  where  several  sources  were  reporting  the 
same  event,  our  source  was  selected  according  to  the  following  priority  list: 
ISM,  PDE,  SDAC/  LASA,  NORSAR. 

The  "comment"  field  provides  information  relative  to  our  pro¬ 
cessing  of  each  event,  and  the  following  codes  have  been  used: 

ND  :  Not  detected  on  NORSAR  data  by  the  TI  analyst. 

MBN  :  Body-wave  magnitude  has  been  computed  based  on  NORSAR  data. 

E  :  The  event  is  a  presumed  explosion. 
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"I  T  Ai 

S IB/014/O’N 

01/14/7? 

C  3  .  /  C  .  ?  C 

67.  CA 

1 71  .r p 

atr 

*4  O 

.  • 

p 

IP A/014/22N 

01/14/7'' 

??. ir.04 

77.^6 

44.fr 

a  n  f 

5.1 

p 

K'JR/Ol  5/00N 

01/15/7? 

go. PP. ?? 

4  o  .  4  A 

1  fc.0L 

MIC 

7.0 

S 

FRS/015/1«N 

01  /  15/7? 

1  *3 . 0  7  .  6  c 

c  7 . 4  A 

1  2  0 . 7  p 

1  3 

4.7 

P 

SI 1/015/20M 

01/15/7? 

?  c .  ?  l . c  r 

4  C .  3  A 

70.  OP 

ADP 

4.4 

c 

SI  2/015/? ON 

01/lc/7? 

20.4  c. 2? 

?0. 7  A 

7?.  nr 

NCR 

4  .6 

D 

K  A w/ 0 1 6 / O^N 

01/16/7? 

04. ?P. 1 6 

4  «:  .  t  A 

16  ?.cr 

MR 

7.5 

5 

K Aw/016/  1 1 N 

01/16/7? 

1 1  .rr. 40 

5  5.65 

16  :.?r 

7  c 

7.0 

<; 

IP  A/01R/21N 

C1/1F/7? 

?1  .  1 ?.  0? 

77. 4A 

4  «.  7A 

A  p  ( 

4.0 

P 

IT  A/0  18/? IN 

01 /IP/7? 

?  7 .  1  6  . 1 7 

44. ?A 

? .  r  p 

2  5 

4.1 

P 

rnn/0?0/0?N 

Cl /20/7? 

07.1^.07 

3  6.65 

?7.1  P 

ATI' 

4  .« 

S 

KUR/022/01N 

01 6??/7? 

Cl .41 . ?4 

50.  OA 

1 C?.0P 

Ant 

4 . 7 

n 

TI1R/02  ?/ 1  7N 

01 /??/?? 

1  7.  17.7] 

77.  4f 

K.cr 

1  7 

4.4 

r 

K  A M/O?  5 / 1 CN 

01 /2C/ 7? 

l  o .  r  ? .  4  o 

c3.o\ 

v  c.°r 

ATR 

4.6 

p 

IT  A/O25/20N 

o  i  /  ?  *  /  7  r 

20.74. 

43.  ^ 

1  «.4( 

7  6 

4.4 

r 

CRF/0?6/l?N 

01/26/7? 

1 ?. P4. 

7  4 . 5  A 

"’5.51 

AT  r; 

4  .0 

s 

K  AM / 0 ?  7 / 20N 

01/?7/7? 

2  0  •  7  7  .  ?  r 

5  5 . 7  N 

16  7.  ?r 

4  C 

3.° 

s 

FC.  S/02P/04N 

0 1  /  ?  P  /  7  2 

C4 .22. 20 

27. c5 

1  ?  4 . 5  r 

Arc 

4.4 

s 

PAK/02P/1 ON 

01/?p/7? 

1 0.26. r4 

?6.65 

4  4.3  4- 

ATI  P 

4  .0 

A 

K IP/0?R/?0N 

01 /?P/7? 

?0. ?f  .  1  0 

47.  ON 

7  7. 0  r 

A  pi> 

4  .6. 

P 

CPS/07R/21M 

01 /?°/7? 

? l . Rc.rc 

44. 0  A 

1  6.  or 

A  UP 

4.0 

t 

KUR/02P/23N 

01/25/7? 

?7.P?. 51 

iti.n 

1 R7. ?p 

NP  R 

7  .r. 

6 

IP  A/Opc/OSN 

oi  /ro/?? 

oo.  l-c.  50 

7  0.05 

4?. OF 

ATP 

7  .P 

A 

KAM/032/ ION 

07/01/7? 

l r. u  .cc 

5  4  .  0  fv 

14?. Of 

Arc 

4.1 

5 

KAM/033/OAM 

02/02/7? 

04.26.50 

5  c  .  75 

16?.0P 

\nr 

7.7 

c 

KIJP/033/00N 

02/0?/ 7? 

0  5  .  r  F  .  c  1 

4  6  .  0  A 

1 6.  4 . 4  r 

A  P  R 

3.6 

r 

) 

KlJ»/03  V17N 

0?/0?/ 7? 

1  7  .  L  6  .  ?  o 

T.  7f 

1 1  0. 1  r 

API 

7.6- 

5 

YIJN/034/07N 

0  ?/()?/  7  ? 

07.??.4O 

?  7 . 4  5 

1 0?.4C 

ATI 

4.4 

P 

IT  A/035/02N 

02/04/7? 

02.4?.  10 

4).  7  6 

1  7 . 7  r 

2  6 

4  .  5 

0 

PA  I  /035/Q3N 

02/04/7? 

03.74. 54 

r  1 . 4  5 

1  1  P.  O'" 

MR 

6»  •  2 

5 

IT A/03S/04N 

0?  /() 4  HI 

C4.4C.r'r 

4?  .  OA 

;  1 ,  ?  r 

A'p  P 

4.7 

P 

IT  A/035/0ON 

07/04/77 

CP. I P. 77 

47  .  OA 

1  ^  .  ?  f 

5  7 

4.4 

P 

IT  A/035/1 7N 

02/04/?2 

1  7.  1  5? 

4  3. 4 A 

1  ?  .  7  r 

?  -> 

4 .4 

P 

ITA/03S/19N 

02/04/72 

10.02.56 

44  ^ 

17.75 

A  r  ;j 

4  .  P 

P 

IT  A/036/01N 

02/C4/7? 

C  1  .  ?  6  .  ?  3 

4  ’I  .  A  A 

1  0 . 7  r 

\ro 

4  .P 

P 

IT  A/036/0  3N 

02 /OP/ 7? 

03.40.4^ 

4  7  ,  ?  A 

1  '  •  7  r 

Apr 

4.4 

P 

IT  A/036/05A 

0 ?/0c/7? 

C  5  .  C  c  .  R  1 

4  7  .  TA 

1  7 . 6  r 

APR 

4.6 

P 

IT A/03A/07N 

0?/0c/7? 

C  7  .  C  P  .  1  ? 

47. 4A 

r;. -r 

APR 

4.7 

p 

AFR/037/01N 

0?/0t  /  7  ? 

01  .74.?? 

4  4 . 0  6 

17.75 

a  r  p 

4 .0 

c 

KA//077/QPN 

0?/06/7? 

C  H  .  C  ?  .  4  7 

44.  ON 

>o.or 

NCR 

4.7 

A 

IT  A/037/21N 

0?/06/7? 

21 .44. ?r 

4  3  .  H  A 

1  7.  ?r 

APR 

4.4 

n 

ITA/0’0/1?m 

0?/0P /7? 

1  2.1  <-.  14 

4  3.75 

l'.-r 

»  r  R 

4  .6 

p 

FKZ/041/0SN 

0?/  10/  /? 

05 .  rr,  r 7 

40.05 

7P.  <  ‘ 

n 

r.  r 

• 

c 

IRA/041/0  ON 

0?/10/7? 

rr .04. cc 

7  .  6  A 

ci'.rr 

AMR 

7 .0 

P 

IP  A/041/1  AN 

0  ?  / 1  C  /  7  ? 

1  C.  1  6 

ir. 

rC.  ^5 

45 

4 . 1 

p 

A-7 


r  iM'VMT 

•J  a  \> 

MR1 

vf 


*  PM 


Vp'l 


Y  nv 


f  r 


TABLE  A  - 1 


E  VFNT 

08  9  IGNAT IHM 

S IN/042/05N 
T IP/042/ 1 2 N 
KAM/0A2/21N 
KUR/048/05N 
GP.  r  /  044/  1  3N 
KOM/044/ 2  2N 
KIJR  /046  /  1 6N 
GP  F  /  04  7  /  00N 
S IN/047/23N 
KUR/049/18N 
SIN/051/1  ON 
PK  H/  0  5  1  /  1  ON 
KAV/Q51/20N 
K AM/p 52/22N 
YIIG/O  5  2 /  2  TN 
MriN/05^/0  IN 
HIN/053/0PN 
KIJR/054/03N 
KAV/054/ ! ON 
KUR/055/1 ON 
KllR/0^5/  1  «N 
KIIR/056/1ON 
WRS/0C6/22N 
KIIR/056/22N 
KUR/057/0SN 
KA^/057/CON 
F&  S/057/ 15N 
VUN/057/ 1PN 
l 0 1 /O  58/0  8N 
LP1/058/10N 
LC1W/058/11N 
L  OM/058/ 1 7N 
PA  I/05P/22N 
KUP / 050/0  IN 
PAK/0C<3/05N 
KAm/ck5/11N 
Arr/059/  1 PN 
*AM/059/2CM 
I PQ/OfcO/OPN 
TP  A/062/1  4N 
Al M/067/1 9N 
K  AM/063/00N 
NS  I/06  3/05N 
KrM/06^/CPN 
VlJG/06  V  2  IN 
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PATC 

no  ir  u 

T  I  MP 

LAT 

1  on 

8IFPTM  8p 

S^URfC 

RLTN 

rnMurr  t 

02/n  n-> 

05. 5r.A6 

39.  9N 

77. 4E 

23 

4.9 

p 

o  ?/u  n? 

1  7  .  ?  r .  4  7 

2  5.  ON 

P  7.  OF 

K'n  p 

/  .3 

M 

02/11/72 

?1 . 36.  17 

56.  1  N 

16?. or 

44 

4.6 

P 

02/13/72 

05. ?4. 57 

43. 5N 

1' 7. OF 

r  pp 

T  .8 

5 

02/13/72 

1  3.07.11 

37.  1  N 

24. OF 

27 

4.5 

P 

02/1  V7? 

?7.36.54 

5  ^ .  2  N 

165.58 

NOR 

0 . 0 

5 

0  2 / 1 r / 72 

16.45. 7? 

45.  ON 

157. OF 

NPR 

8.1 

4.5 

5 

02/16/72 

00.4?. ?4 

36. 9N 

?8.nr 

N  np 

p 

0  2/) A / 7  2 

?  ^  .  1  Q.  ?C 

41 . 7E 

FO.  79 

2  0 

4  .  r 

p 

02/18/72 

15.02.34 

63. 6N 

187. HF 

36 

4.7 

S 

I 

02/20/7? 

10.22.46 

3H.5N 

00. 5F 

16 

3.9 

VPN 

02/20/7’ 

10. OP. 46 

47.  ON 

148. Of 

07 

4,7 

1 

02/2C/72 

70. ^6 . 1 1 

50.  «N 

18] ,8F 

NPR 

4 . 1 

4 . 5 

1 

I 

02/21/72 

7?.cc. 

54. 4 N 

161.38 

NPR 

t 

02/21/7? 

? 7.  C2.  c5 

41  .ON 

??. Ti¬ 

NPR 

8 .0 

I 

02/22/72 

Cl  .  5?.  76 

4  9.  ON 

ll  5.08 

NTR 

8 . 1 

I 

02/22/7? 

OH. 14. 26 

36. 6N 

6  8 . 8  f 

N'CP 

4.0 

I 

M  PM 

02/23/7? 

07.42.41 

43.  9 N 

14  8.3F 

7  n 

4 .0 

T 

07/23/77 

l g .37. 29 

55.  ON 

163.08 

NPP 

3.7 

I 

NT 

02/24/7? 

1 0. 1 c. 77 

4  H  .  .3  N 

155. 7F 

NPP 

5.0 

I 

02/24/7? 

16.17.  36 

49.  ON 

15  P.08 

NCR 

T  .5 

I 

NO 

02/25/77 

1G.5C.  ?<5 

46. ON 

187.08 

NPR 

3  .  « 

I 

02/28/72 

22.34.49 

50. OK 

TP. OF 

NPP 

3 . 7 

I 

02/25/72 

77.47.C7 

49.  ?N 

1 56. OF 

Nor 

4 .0 

I 

0.7/26/7? 

C5.«=P.  2? 

48.  p N 

157.68 

NPR 

4 . 0 

I 

02/26/72 

0C.C6.7? 

55. ON 

167.08 

NPR 

3 . 3 

I 

T 

NT 

M  PN 

02/26/77 

1  6. C6.  *  2 

53. 3N 

13P.7F 

NPP 

3  . « 

T 

02/26/7? 

l  p.  56.13 

27.  1  N 

100. op 

NCR 

4 .7 

l 

02/27/7? 

CP. 42.  5 c 

PH.  ON 

74.  OW 

NOP 

3.3 

I 

NP 

02/27/72 

10.07.C3 

87.  ON 

63. 5h 

NCR 

8.Q 

I 

02/27/77 

1  1 .07.  l  9 

so.  or 

95. Oh 

NPR 

?  .  5 

I 

0  2/ ?7/ 72 

1  7.*--C.  25 

P6.2N 

77. ?h 

NPP 

4.4 

I 

02/27/7? 

??. 1 5 . C  3 

55. ON 

03.28 

NPP 

4 . 5 

I 

02/78/7? 

01 .04. 2? 

46.  ON 

1 4  P,  08 

NnR 

4  .  ? 

I 

NT 

02/28/7? 

Oc. 1 P. 56 

3  6.  7N 

‘  1  .4F 

Nnp 

4 . 2 

I 

M  PM 

02/ 2  P/7? 

11.35. 31 

58.  0  K 

163. OF 

NnR 

4 . 1 

I 

NT 

02/26/7? 

1 ^ 12. 75 

3  4.  ON 

6P.7F- 

NPP 

4.4 

I 

WPM 

02/28/7 ? 

20.C6. CO 

5  5 . 1  N 

164.78 

NOR 

3 .4 

I 

Nn 

C2/?«/7? 

C-7.C2.51 

3  2  •  3  N 

88. 6r 

NTR 

4  .0 

I 

PPM 

07/02/7? 

18.1 C. 13 

31  .  6  N 

4  7.18 

NPR 

4.6 

I 

vpm 

03/02/7? 

1  5.57.4? 

43. ON 

76. OF 

Nr  R 

T  .5 

I 

VPN 

07/03/72 

P0.?r.2? 

57.  p\ 

15  0.78 

N°R 

8 . 1 

1  .  R 

I 

07/07/7? 

C  5  .  :  f  .  9  3 

77.5  N 

1  !  4. 78 

NPR 

I 

03/07/7? 

OR. 1 3. cc 

6  5  .  P  N 

163.98 

NCR 

6 . 1 

I 

C3/03/72 

7  1  .  2  T  .  8  1 

8  4. 7N 

18.48 

7? 

4.0 

I 

A-8 
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EV  ENT 

OP  I  G  I  N 

DESIGNATION 

CATE 

TIME 

KUR/06 ’/23N 

03/03/7? 

2  VIC. 41 

S1N/064/04N 

03/04/7? 

C4.00. CS 

KAM/066/06N 

o?/5fc/7? 

C6. C5. CP 

KUP/066/CSN 

03/06/7? 

0  9 . 5  «  .  C  0 

DKH/066/ 1 9N 

0B/C6/72 

15.1  3.  24 

CHI/066/23N 

0’ / C  6 / 7  ? 

??.  1 7.5’ 

YUG/067/05N 

03/07/7? 

08. ?1 .21 

OK  H/06  R/02N 

O’/O  P/72 

02.  ’P.l  1 

IR A/068/21N 

03  /OP / 7? 

?1.4C.ii 

PUL/O60/22N 

O3/0P/7? 

22.04.0? 

KA7/070/04N 

03/10/7? 

04.56. 57 

KUR/07C/06N 

03/10/7? 

C6. cc. 1 P 

APC/071 /06N 

0  VI  1/72 

06.47. 07 

KAS/071/13N 

0  V  1 1  /  7  ? 

n.?i.  ’r 

KU°/073/0?N 

0 1  /  1  3  /  7  ? 

0?.l 1 .cc 

AFC/073/05N 

OB/13/7? 

05. 4 c.  i  ? 

T  1  P/073/  1  0N 

03/1V7? 

1  P.  ?7.  07 

TIR/075/06N 

0’/lr/7? 

06. CC. B ’ 

KUP/C77/07N 

OB/17/7? 

07.if.C? 

T  AD/077 /0QN 

03/17/7? 

05.17.  ’  1 

IP  A/077/  1 7N 

03/17/7? 

17.11.?a 

KAS/077/73N 

OB/17/7? 

2’  •  3  "V 37 

KA7 /07P/C7N 

0° /l P/7? 

07.11. 5C 

KAM/070/1 3N 

0’/l«/7? 

1 ’.r2.  14 

KAM/O70/1PN 

0  .3/1  »/ 7? 

1  P  .  2  5- .  ■»  7 

OK  h  /0 7  P  /  1  6N 

0  3/1  P/72 

15.1 7.2C 

K  IR/O70/19N 

O^/lP/72 

15.  c  4 .  ]  P 

CMI/079/03N 

0’/  19/72 

C’.34. 31 

KUR/0«0/lAN 

0B/20/7? 

14.  CP.  1"* 

F l N/OBO/ 2  IN 

03/’O/7? 

?1  .47.^ 

CK7/1 53/01N 

06  / C 1 / 7  2 

01  .? s. ?6 

MON/153/11N 

06/01/7? 

1  l.??.lc 

GP  F/153/ 1  BN 

06/01/7? 

13.44.11 

KP8/153/21N 

04/01/7? 

21  .4?.45 

IR  A/154/00N 

06/02/7? 

00. 1 ?.  1  ’ 

KIJR/154/01N 

04/02/7? 

01  .  c  3. C7 

SI 1/164/ 04N 

06/07/7? 

04.21 . 45 

S  I  ?/  1  6  4 / 0  4N 

04/02/7? 

04. ??.  16 

SIN/154/0 5N 

04/02/7? 

05.  11  .  1  -> 

SIN/154/06N 

Q6/0V7? 

0  V  ’  C  . 

TSI/1 54/16N 

06  /02/7? 

1 ' .4C. ?? 

CHI/184/20N 

06  /O?  /7? 

?  0  .  ?  ?  .  5  c 

RYU/  1  6«=/0?N 

16/03/7? 

0?.  1  '  .  51 

IP A/1C5/0FN 

06/0' /7? 

o  •> . 0 1 . ‘  r 

IP  A  /  1  c 5/0  BN 

06/04/7? 

0 V ’7.45 

l  A  T 

i  r  n 

nr  pth 

MP 

S  ni|PO  c 
R|.  T  M 

COMMENT 

50.25 

1 5c.7E 

NOR 

4  .  c 

I 

4  0  •  ?  A 

7  5.05 

NOR 

4.6 

I 

63.  5N 

1  6  0.  5  r 

N'5R 

3.9 

1 

r 

45.  OK 

1  C0.0F 

N  r  p 

7.7 

I 

56.  ON 

140. OP 

NCR 

4.? 

1 

40.  ON 

103.05 

NGP 

4.6 

I 

NO 

42. ON 

?1 .05 

M’P 

?  .7 

I 

►  i 

M  .  ?N 

1 f 1 .or 

NJ-'P 

4.? 

I 

?  7 . 6  N 

c  6 . 7  c 

46 

4  .O 

I 

40.  PN 

?  ?  .  Q  F 

Mf  P 

’  .6 

1 

M  0\j 

45.  PN’ 

7  p  .  ?  r 

0 

6.4 

1 

r 

4  c  .  1  N 

1  4  5 . 6  c 

NCR 

3.7 

I 

HI  .  OF 

1  r7.0r 

Nfi  P 

4.7 

I 

Nr 

•a  c  .  OM 

76.  OF 

N  IP 

4.1 

1 

V  U  M 

45.  ON 

1  c  3  •  Of 

\ro 

’  .  P 

1 

37.  ON 

70.  OF 

NOR 

4  .0 

I 

•*  MV 

34.  Or 

r'  3 . 0  F 

N  r  p 

4  .  1 

1 

V  PM 

3  0.4  N 

F4.6F 

Nrp 

r,  i 

• 

1 

>  <  P  ►  | 

4  Q  .  0  A 

1  5  £.  .  ?  C 

n  r  p 

6.2 

I 

43.  1*1 

6  ° . 7  r 

’6 

6.? 

1 

?  VON 

c  4  . 

N<  n  r 

*3  rs 

• 

1 

.'PN 

32.0N 

’ c .  OF 

NOP 

7.6 

T 

*  [1! 

4  7.  ON 

PI  .or 

n  r  D 

3.6 

1 

MOM 

*•7.  ON 

1  ‘-3.0r 

NPP 

’  .6 

1 

5  0 . 6  N 

1  5  6.  7  r 

\r  i> 

4  .7 

T 

54.  ON 

150. OF 

Nf  t 

’  .7 

I 

41  .ON 

7?. OF 

NOP 

’.? 

I 

'VN 

4?.7fv 

’B.  1  r 

NFR 

7.5 

1 

47. (}N 

7 64. OF 

4  .0 

1 

'Ml’  1 

40. ON 

»  o .  o  r 

N  r'  P 

7  .<> 

1 

7,  P  *  1 

5?.  ON 

70.  05 

Nf 

*3  t , 

• 

S 

v  r 

44.  ON 

l 03. or 

vrc 

“  •  7 

N 

75. ON 

?4.or 

N;r  p 

4.1 

5 

55.  ON 

1' 4.0r 

fT« 

’.7 

c 

35.  ON 

5  7.  ()r 

N  r  c 

6  .1 

N 

6  0.  ON 

] r?.0r 

m  ■”  p 

-J  3 

• 

5 

'  r 

4  ?  .  G  N 

p:.or 

N  r‘  P 

’  .f 

K* 

4  ?  .  0  f 

j  >.  Or 

N  r  c 

3.7 

4 

4  ’  .  0  f 

■’1  •  5  f 

N  r  P 

7.6 

N 

4  ?  .  0  * 

’  1 .0r 

Mr  p 

*7  r> 

• 

N1 

7  6 . 0  N 

r  ?  ,Or 

\r  r 

M 

N 

?  r .  4  m 

5  6  .  Q  f 

Np  f 

4  .  ’ 

0 

2  3 . 5  N 

1  ?  r  .  c  u 

N  r  i 

r,  ^ 

• 

r> 

?  5 .  n  n 

r  %  .  Of 

\r  i.' 

• 

*. 

’0.  ON 

‘  4.  V- 

^  '  j 

‘.7 

N 

FVpNT 

f'f  S  IGN^T  IflN 

KAM/l«;t/07M 
KPM/156/12N 
KPM/156/l 3M 
THP/166/1  6N 
TS  1/156/23N 
K AM/IS7/04N 
GPF/1  6 7/1  ON 
IP A/l*7/llM 
PAK/ 1=7/ 1 1M 
Nr  J/15  7/nN 
K»I»/15P/0?N 
•UIR/1  6P/06N 
K AM/163/1  ON 
F*  /  /  1 Sc /O IN 
Kl>*/1  59/06N 
TA  1/ JfcO/OOW 
K  A/16C/C9N 
T Al/Jf O/ION 
THR/16C/12N 
PUP/1 60/16N 
TA  I/ltO/lfcN 
C  AM/ l 60  / 1 7N 
T  IP/HC/P’N 
KUP/161/OON 
ru  F/1M/07M 
Tl|P  /  1  6  1  /  l^M 
IPA/163/03N 
PAK/UP/UN 

KAV/1(,?/14H 

<IJ1/16?/?3N 
KU?/l/>7/?3N 
KIIP/1FA/OON 
IPl/l‘A/1 3N 
|o?/164/13N 
1° A/ 165/0 ON 
* AM/i 6S/04N 
r  AW166/00N 
I  °  A/166/046' 
IF  1/166/1 2N 
IP  7/166/1 ?N 
I  T  A  /  1  6  fc  / 1  PN 
GP  r/ 1 67/00N 
*0*/  167/  1?N 
T'|0/167/14N 
Ki^/ ) 6°/ C9N 
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">  ATF 

ni>  |  o  i  * 

TTMr 

1  AT 

1  rN 

op  PTH 

pn 

cnii'CF 
«»|  TV 

roMMfvT 

06/04/7? 

01.6?. ?o 

6  3.06 

1 60. OF 

6  V  P 

4  .0 

S 

06/Q4/7? 

1  ?  •  *•  7.  3  0 

63.06 

160. OF 

NOP 

7.4 

5 

o*>/n^  ni 

» J.0?.C7 

6 A.  06' 

1>  5.  Or 

NOP 

?  .6 

c 

06  /C4/7? 

16.?4.34 

36.^6 

:a.?f 

17 

4.1 

p 

06/06/7? 

??.??. 1 P 

27.06 

97. OF 

NOP 

3.6 

N 

1 

06/05/7? 

06. 1?. 56 

56.  ?6 

1  A  3  .  1  F 

NOP 

6.3 

P 

06/06/7? 

10.4A.6C 

37.  AN 

?1 ,4F 

60 

4.? 

0 

r 

06/05/7? 

1 1 . 1 7.67 

7 A.  f)6 

46.  OF 

NOP 

7.‘? 

N 

06/0K/7;> 

1 !  .  c?.  *•? 

?S.  66 

7C.3F 

?  7 

4.0 

P 

Li 

06/06/7? 

19. OC. 1 ? 

P6.  66 

?  fi  .  5  F 

Nr.” 

A. 6 

5 

06 /C6 / 7  ? 

n?.0A.44 

4A.06. 

14  3. 0  r 

6C° 

’.A 

6 

n 

06/06/7? 

06.32.10 

40.  06 

165. OF 

NCR 

7.6 

6 

\r  L 

J6/C.4/7? 

10.4?. 3? 

6  =  .^6 

1M.3C 

NOP 

4.7 

P 

06/07/7? 

nl  •  ?  7  .  r-  7 

4  5  .  fl  6' 

79.  ?r 

0 

».  .c 

P 

r 

06  /07/7? 

of. cc. ?o 

56.06 

166. OF 

6*0  P 

0.0 

6 

Li 

06/00/7? 

G=>.1  A. OP 

21.16 

l?0.?r 

6|OD 

6.4 

P 

06/06/7? 

0  c  .  0  r  .  ?  1 

74.16 

46.  ?r 

1  6 

4.C 

P 

ri 

06 /0° / 7  ? 

1  0.  1  7. 4  A 

71.06 

\,?0.?  = 

NOP 

4.9 

D 

06/0P/7? 

1  <1.  A6.  1  = 

41.06 

44.  oK 

6TP 

4.1 

S 

Li 

06 /OP/ 7? 

16.CP.C6 

15.06 

=  4.0  = 

•ynp 

4.7 

N 

06/03/7? 

1  =■  .  A  4  .  ?  4 

?  1  •  06 

1  ?  0. ?  F 

MOD 

4.7 

0 

r  ] 

06/0°/7? 

17.75.5? 

43. ?6 

A  7  .  ?  F 

NON 

4.6 

r 

LI 

06/06/7? 

??. 1C. 1? 

?°.  66 

4,?.  3  = 

44 

4.7 

0 

06/0c/7? 

00. 16. A? 

4  7.06 

1 =?. Or 

Noo 

4  .4 

s 

06/05/7? 

07.4?. ?0 

74.76 

?►.  =  = 

6T  0 

4.0 

p 

0* /0«/7? 

1  =.A?.?7 

77.06 

a  4.  or 

6jn  p 

4.0 

N 

C6/10/7? 

0? .  ?5. ?? 

71.06 

rl.n= 

6trip 

3.6 

6 

r~7 

C6/10/7? 

11.75.11 

?  M.  ?  6 

►  4.6  = 

6)0  Q 

4.5 

P 

06/11/7? 

16. 1A.01 

63.06 

160.0= 

'IMP 

?.7 

6 

NT 

C6/11/7? 

??.??. 04 

40. 06 

16?. 0= 

NOP 

4.0 

c 

06/11/ r; 

?».??. 44 

4  7.06 

Is?. OF 

6’”  3 

4.7 

5 

O^/l?/7? 

CC.lc.16 

44.06 

1' 9.0= 

60  U 

7 .7 

5 

0  6  /  1  ?  /  7  ? 

l?. ’4. Cl 

3  7.16 

4f  .  ?  = 

NOP 

6.4 

D 

06/1?/?? 

1  1,U/C 

7  7.16' 

46. ?r 

N  ?  P 

6.1 

P 

pj 

06/17/7/ 

0  C • 6  c .  37 

77.16 

4  ►  .  ?  F 

?7 

c.l 

P 

L- 

06/1 ?/7? 

04.63. 7C 

6C.06 

16?. 0= 

NfP 

3.o 

5 

06/16/7? 

00.4=.  =  A 

40.16 

=  1  ,oc 

47 

4.7 

n 

06/16/7? 

04. ? A. ?o 

7?.  06 

44.  1  = 

6-nn 

6.? 

p 

06/14/7? 

l?.ll.?3 

71  .C6 

r  r .  0= 

MPR 

?.6 

F 

LJ 

06/14/7? 

1  ?.  ?•-.  C  = 

?  7 . 0  6 

c  6 .  OF 

NOP 

?.6 

\ 

06/14/7? 

1  °.  6*-.  6? 

4  >.  76 

1  7.4  = 

14 

4  .4 

P 

r 

06/16/7? 

0  ).?’.?4 

?  V  .  ?  6; 

??.?  = 

?6 

4.0 

p 

06/1=/?? 

1 ’.49.1 J 

=  4.^6 

' * p.OF 

6flO 

?.6 

5 

06  /  lr/7? 

14.1C.0? 

3  =  .  0  6 

?  =  .0  = 

»nu 

7.3 

• 

r 

06/1//7? 

C5.a4.ai 

=  4.  C6 

1 '  1  .or 

V"»D 

A.l 

6 

*-10 
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FVFNT 

0F$  IGNAT inM 

HIN/160/16N 
KAM/160/72N 
IKA/160/23N 
AUS/169/09N 
KIJR/169/19N 
TIR/17C/04N 
KUP/17C/05N 
Kli;’/I71/I«N 
KUP/171/22N 
IPA/172/05N 
FPft/  172/09M 
KAS/172/15N 
KAM/173/0CN 
KM  2/ 1 73/OCN 
TUP/173/05N 
KA«*/173/lON 
TUP  /  I  7  3/  1  AN 
ITA/173/15N 
KUR/174/02N 
TUP  / I 75/04N 
GRF/175/07N 
T  AT/ 1 75/  16N 
IPA/176/C6N 
VUG/1 76/07N 
MIN/ 176/16N 
TAD/17A/1RN 
VUG/177/04N 
HIN/177/07N 
KAM/177/17N 
TA  I/178/0RN 
K  A**/  1 78/ 1 7N 
MIN/178/2  ON 
AFG/179/05N 
P  AK / 1 79/06N 
KAM/17R/06N 
9IIP/179/09N 
WP  S/ 175/ 1 2N 
MIM/179/15N 
YIIG/1AO/01N 
TSI/180/03N 
KPM/1 R0/06N 
CYP/lflO/ORN 
II AP  /  1 P0/05N 
KAM/1R0/14N 
CKZ /  181/00N 
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r»ATc 

no  1  r.  j  n 

TIMf 

LAT 

LTK 

r.c  pTH 

»<« 

ft  PI  l po  C 

»IT» 

rr.v  1 

06/16/7? 

IP. 67. ft? 

36.  "K 

65.56 

40 

4 .4 

P 

Oft / 1 6 / 7? 

22. 12.15 

ft  5 . 0  ft’ 

1*7. or 

K!°P 

5  a 

. 

5 

vr 

06/16/7? 

21.5?.?7 

34.06 

46.06 

KPt 

3.7 

K 

06/17/ 7? 

C6.C?. 4P 

4«.n 

1 4  .  ft  r 

KC  5 

4  .A 

D 

06/17/75 

1 Q. 1 °.  21 

44.  ?W 

lAO.jr 

6  4 

4.6 

P 

06/1 P / 7  7 

04.20.47 

33.06 

C3.pt 

k  ro 

4.5 

K 

06/1 P/7? 

CP. 1 C. c4 

4".  CP 

1  ft'-.or 

«g(m 

5  .O 

ft 

ft/l°/52 

!  0  •  C  7 .  *■  5 

43.  *K 

1M  .ftc 

Kuo 

4. ft 

p 

06/19/73 

5P./-1.4’ 

40.  OK 

1ft  7.. Or 

vr.R 

4.1 

s 

ft  p 

C6/20/7? 

Cc • 1 7. 45 

79.0V 

fti.or 

\rc 

4.0 

Kf 

06/50/72 

C5.1P.CC 

ft  ?  .  C  V 

151  .OP 

K  0  K 

3.7 

ft 

KO 

06/20/75 

15.34.57 

5  5.  OK 

7ft.  0  = 

NCR 

5  t 

• 

ft 

o*>/?i/y? 

CC . 1 2. ft  B 

ft’.  OK 

161 .0C 

\i"r 

*  •  " 

K 

Oft/?i  n? 

OC. 15.05 

F4.0K 

1  so.nc 

3.5 

ft! 

Oft  /  ’  1  /  “’2 

C6.C6. 17 

4  0. /ft 

5  0.  Or 

•cnc 

<  .1 

P 

06/21/72 

10. 4?. 4ft 

54.  OK 

161. OF 

K  TP 

4.3 

ft 

06/21/72 

I4.r5.rc 

57.  Oft 

41  .or 

\rp 

3.3 

ft 

Oft/21/7? 

l  c  •  06  .  *>3 

43.  “K 

1  3.5- 

4 

4.4 

P 

0ft/2?/75 

0?.?ft.ftl 

40.  OK 

1  44.0c 

vnu 

4. ft 

N 

06/52/7? 

C4 . 2  ft. ?  7 

41.06 

5  0.  or 

KTO 

-  1 

• 

ft 

06/27/7? 

C5. IP. tft 

57.  ?fti 

5  1  .or 

ur'p 

’  .4 

ft' 

06/55/7? 

16. ftc.4P 

37. CK 

7ft. or 

vr  n 

3.7 

K 

06 /2ft / 72 

C6.ft7.C7 

23. CK 

ftA.or 

ft  flp 

3.  ft 

ft 

06 /2ft / 7? 

07.17.ft6 

43. 7K 

1  6.rr 

ftoj. 

c  7 

• 

P 

Oft/24/72 

16.14. F4 

36.  OK 

6  0.  pC 

M'.r 

3  .4 

ft' 

06/24/7? 

1 1 C 

55.  ok 

74.  V- 

K;? 

3.4 

ft 

06 /?  e / 7? 

04.55. 1' 

44.  OK 

1  ft,  dp 

ftr'P 

4.4 

P 

06/5R/7? 

07.  f  ft.ftft 

3  4.7ft 

65. 6r 

4  A 

4.7 

D 

06/2^/72 

1  7.5ft. ftQ 

ft  4 . 0  V 

160. Or 

KfO 

4.1 

ft 

KT 

06/56/72 

0 1 .  r  p .  ?  ft 

21 . 1  ft 

1 ?0. 5C 

KTO 

ft.o 

P 

06/26/7? 

1 7.?:. 57 

ft 6  .  Oft.' 

10 J>0r 

upr 

3.6 

s 

ft  r 

06/26/72 

20.  ftC.ro 

3  6.  QK 

fto. or 

KPP 

3.7 

K 

Oft/57/7? 

Cft.c7.45 

3  °.  OK 

6ft.  Or 

K  O  C 

4.0 

K 

06/27/7? 

C6.5C.«4 

?  r  .  7W 

50. 5  C 

l? 

•*’  .ft 

P 

06/27/7? 

C6.<,f.C5 

ft  4 . 0  K 

l=o.or 

M'? 

3.=> 

ft 

Oft/. 7  7/7? 

c**.  c  ft .  ft  ? 

56.  ?, 

06.6- 

•> 

4  .4 

f> 

Oft/27/ 7 5 

12. 50. 5b 

c!  .  Oft 

6  7.  )C 

k  re¬ 

3.6 

ft 

0*>/?7/72 

1  ft. cc. 3e 

36.  ^ft 

AO.  ftp 

ft  3 

r.l 

P 

06/2P/7? 

Ol.45.ftA 

*  3 . 0  K 

5  0.  HF 

\r  0 

4.C 

!) 

06 /?  P / 7  5 

C3.Cr.ft5 

51,0ft 

r  1  .  Oc 

KPR 

*  .6 

ft' 

Oft /?°/ 75 

C6.CC. 2? 

ftft.Oft 

164. nr 

ftr p 

1  .4 

ft 

K  o 

06/23/7? 

CP. 1 6 • ftc 

3ft. Oft 

3?.  nr 

K('P 

4.3 

K 

06/2F/75 

C5.4c.?ft 

?  7.  *  K 

5  3.  «r 

1  ft 

ft.  4 

F 

Oft/?o/72 

14.  ftf. 

ft  3 . 0  K> 

161 .Or 

•  1  r**  r 

;  o' 

•  ‘ 

ft 

06/20/3? 

no.4i . n? 

ft4.ov 

t ■  7. or 

1,7 

ft 

A-l  1 


TABLE  A  -  ] 


rVFMT 

0 F  $  IGNAT  f  (IN 

A  FG/ 1 R 1 /O^N 
lfWlR;>/17N 
T  A  I  /  1R  2/  1  RN 
I  R  A  /  1  R  2  /  ?  ON 

TR  A/1R4/12N 
T P A/ 1 «4/ 14N 
I  0  A/ 1 R6/02N 
AFG/  1R5/03M 
l» A/1HF/17M 
I D A/185/1 ?N 
IP  A/1PS/71N 
K 1 1 P / 1 R6/04N 
T *  I  p  /  1  P,'./OftM 
IR  A  /  1  R6/09M 
K AW/ 1«6/ 1 3N 
KIJB/  lfft/  ?1M 
I  R  A/1  87/01N 
SIN/1P7/01M 
5IN/1F 7/07M 
S I N/ 1 87/0 AN 
IRA/1 87/pn»j 
I R A/1R7/16N 
GRF/1R7/1PN 
FK7/18R/01N 
PAK/  188/ 16N 
KtJR/1  RP/1  9N 
K  AM/ ] pc / 0 8N 
PIIR/1  89/]  2M 
R7F/1R9/->3N 
GRF/190/0RN 
KUR/10C/08N 
KIJO/]  Q0/71N 
WRR/1R1/07N 
Mpr/ini/nN 
RYll/lop/O  ON 
R  Y<J/  1 R  ?  /  0  3  N 
R1N/1Q?/10N 
APG/193/0AN 
K(IR/197/0^N 
K AM/ I93/08M 
I  D  A  /  'j  9  ?  /  1  5N 
I°A/19  V22A' 
KIJR/1QA/00N 
P  A  K  /  1 9  A  /  0  1  N 


0«T  F 

06  /  2 c  /  7? 
06/70/77 
06/30/72 
06/10/72 
07/01/72 
07/07/7? 
07/0  V7? 
07/Q3/7? 
07/C3/77 
07/07/72 
07/03/7? 
07  /O’/^ 
07/06/77 
07/04/7? 

07/04/7? 

07/04/7? 

07/06/7? 

07/Oc/7? 

O7/0c/7? 

07/09/7? 

O7/08/7? 
07/06/7? 
07 /06 / 7  ? 
07/C6/72 
07/06/7? 
07/06/7“’ 
07/06/7? 
07/07/7? 
07/07/77 
07/07/7? 
07 /OP/ 7 7 
07/0F/7? 
07/08/7? 
07/09/7-> 
07  /CO/7? 

0  7/  10/7? 
07/10/7? 
07/10/7? 
C7/H/77 
07/11/7? 
07/H/7? 
07/11/7? 
07/11/7? 
07/12/7? 

07/17/7? 
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OR  t  n  T  A' 

T!  vr 

l  AT 

in  a 

OFPTM  pp 

0  •» .  J  ?  .  1 1 

3°. 

71.47 

63 

4.9 

I  7.  *c.  ?-> 

27. 2A 

66.  pr 

Npo 

4.6 

1  Q.  87.  4? 

24. 3A 

121 .  ir 

APR 

4.9 

20. 31 . 3? 

.30.  OA 

63.  OP 

AP& 

4  .0 

02. 10. IP 

54. OA 

166.  OF 

A' OP 

3.4 

12.6f-.07 

30.  1  A 

6  9.Pr 

31 

6.4 

16. 06. C  6 

30.  ON 

6  0.  c‘- 

31 

4  .( 

02.10.00 

30.  1  A 

c0.  P F 

3  f\ 

6.0 

03. R?. 60 

R  6 . 2  N 

71  .  IP 

1  2  P 

4. 3 

1  ^.^1.05 

30.  OA 

6  7.  OF 

NOR 

4  .0 

I*.?*.?? 

32.  OA 

60.  nr 

AOR 

4.0 

21.RF.2? 

RO.  OA 

61  .OF 

4  3 

8.1 

0  A  .  4  ?  .  04 

49.  o A 

1 46. OP 

NOR 

3 . 7 

06.17.  ?r 

41  .  ON 

RR.or 

AOR 

3  .4 

0R.?p. 07 

2  R .  0  A 

6  4 . 0  F 

APR 

3  .° 

1  ‘‘.6?.  ]C 

65.  OA 

1 63. OP 

NOP 

4.4 

21  .47.67 

4  9.06 

1  61  .or 

NOP 

R  .6 

01 . 04.46 

?  3.  ON 

c 4 .  nr 

NPR 

3  .P 

01 .09. 8 3 

4  A  .  <  N 

pi  .  1  r 

A  P  P 

4.6 

0  ? . 4  1.  6  4 

44.  OA' 

°  4 . 0  p 

AOR 

R.6 

O4.0<-./<- 

4  7 . 0  A 

R  7.  or 

APP 

4.3 

09. Rc. 09 

R  R  .  0  r\ 

60.09 

A'OP 

7.4 

16.76.77 

31  .ON 

6  ?  .  0  P 

NOR 

4.0 

1 fi. . CO 

R  7.  OA 

20. OF 

ACR 

4  .4 

Cl. 02. 6p 

4  9.  7AI 

7  R .  0  P 

0 

4.4 

1 6.06 . ’2 

ro.  n 

69. 7f 

C  7 

4  .  R 

1 a.C2. 20 

44 . 0  A 

146. OF 

Nrp 

3.9 

06.17.06 

56  .  OA 

1 63. OP 

AOP 

3.7 

1 7. 06. 12 

2  0.  6  A 

OR.  1  p 

27 

6.0 

2R.43.41 

R?.  ON 

107. or 

AOP 

3.7 

06.46.  1  1 

4?.  OA 

?4.  or 

APP 

4.7 

0  Q  .  2  ^ .  7  7 

46.  1  A 

164.69 

NPR 

4.o 

21  .07.27 

4R .  0  A 

151.09 

A!  PR 

4  .  ? 

07. OC. OP 

82.  OA 

31  .  OP 

NOP 

4.6 

1 R.21.2? 

RF  .  0  A 

1  9.  OP 

NPR 

4.0 

00.6 1 . ?P 

2  9  .  o  N 

131.09 

Nnp 

3.R 

O’. C2. 0? 

RO.  OA 

i?9.or 

AT  R 

3  .  R 

1  c.  03.  3  ^ 
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